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Despite a relatively large list of suspected forms, parthenogenesis has 
been indisputably established in but two species of mites,! namely Tetrany- 
chus telarius L.,? and Pediculoides ventricosus (Newp.).* In view of these 
facts additional cases are of particular importance. An investigation of 
Pediculopsis. graminum (Reut.) seemed of especial interest, for in this 
species Reuter described the process of ‘‘merokinesis,’’* which has found 


no parallel in any other organism. The results of a reinvestigation of 
merokinesis and a study of sex determination in this mite will be pub- 
lished elsewhere. 

Specimens of Pediculopsis graminum (Reut.) were collected from heads 
of Herd’s grass (Agrostis alba L.) at Flushing, Long Island, during the 
summer of 1936. From these collections were established the laboratory 
stocks upon which the present study is based. Not only do these 
mites possess the characters published as peculiar to the species Pediculopsis 
graminum (Reut.), but Enzio Reuter has determined and recorded® speci- 
mens from New York State as identical with the species he described from 
Finnish grasses. For these reasons it is believed that the mite under 
consideration is the same species as that which formed the basis of Reuter’s 
classical investigation of merokinesis. 

Following attachment to the source of food, the hysterosoma of the 
female swells with fluid from one- to five-hundred times the original body 
volume of the mite. Within this taut sac the eggs are formed and develop 
through larval and nymphal stages to adult mites. Mass birth, marked 
by the breakdown of the hysterosoma and consequent release of the sexually 
mature adult mites, is the rule. Occasionally an early larva or nymph 
may be extruded from the mother before such a breakdown; males, espe- 
cially, may gain exit before the mass birth of the brood. It was early 
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found that the unborn mites normally hatch from their embryonic en- 
velopes before the maternal hysterosoma has burst. 

The entire brood dissected from 17 pregnant females totaled 1643 
individuals. Of these, 74 or 4.50% were males. Similarly 11,446 off- 
spring were obtained from 72 pregnant females. Here 511 or 4.46% were 
determined to be males. This confirms Reuter’s observation that the 
males constitute about 5% of the normal brood.® It is interesting to note 
here that every brood contains at least one male. In each of 93 pregnant 
females, no one of which bore more than 10 offspring, there was found to 
be one male while the remainder were females. Nineteen of these mothers 
contained but one embryo, and in each case that one embryo was male. 

To test the mite’s capacity for parthenogenesis, unborn but hatched 
females were dissected from the hysterosomae of a number of mites and 
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FIGURE 1 





a, Metaphase of the primary odcyte division (4410). 6, Somatic metaphase plate 
in cleavage stage of haploid embryo (7050). c, Late prophase of somatic mitosis 
from limb-bud of advanced embryo (7050). d, Somatic metaphase plate in cleavage 
stage of diploid embryo (7050). 


placed in isolation culture. Such isolated females produced young in 
large numbers, a fact which would seem to indicate parthenogenesis. 
The broods of these mites, however, were found to be mixed, the sex ratio 
being close to that observed in broods from non-isolated females. A total 
of 4385 offspring was dissected from 30 of these isolated females. Of 
these, 226 or 5.15% were determined to be males. Each of the isolated 
and dissected females had contributed at least one male to this total. 
The explanation of these findings was found in the mating habits of 
the animals. Reuter has described the copulation of these mites as taking 
place on the host plant.6 Bradbury states that mating often occurs on 
the mother’s abdomen,’ as has been described for Pediculoides ventricosus 
(Newp.).2 In the Flushing stocks of Pediculopsis graminum (Reut.), 
however, copulation may be seen to take place while the mites are still 
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within their mother’s body. The adult mites may easily be observed, 
through the transparent walls of the hysterosoma, to be scrambling about 
and copulating a short time prior to their birth. If females in late stages 
of pregnancy are plunged into boiling water, a pair of offspring fixed in 
copulo may occasionally be dissected from among those massed within the 
hysterosomae. The copulatory position is similar to that figured by 
Brucker for Pediculoides ventricosus (Newp.).® 

Duration of copulation has been observed to vary from 5 to 35 seconds. 
Females have been seen to mate with a number of males, while a male 
has been noted to mate with the same female twice in succession. In 
only four instances has copulation been observed in these cultures outside 
the maternal body, in spite of the fact that males and females are often 
to be found wandering over the abdomen of the pregnant female. In three 
of these four cases, copulation took place within the drop of fluid exuded 
from the collapsed hysterosoma. It has been noted that the mother is 
generally dead when her young are born. It is quite probable that the 
destruction of the tissues of the mother by her active offspring is the cause 
both of her death and of rupture of the hysterosoma. 

Obviously copulation within the mother’s body had prevented isolation 
of virgin females in the first experiment. Accordingly a second series of 
isolation experiments was performed. This time unborn females which 
were still within unruptured embryonic envelopes were isolated. Many 
of these females hatched within the isolation dishes. As Monterosso has 
shown for Pediculoides ventricosus (Newp.), initiation of the swelling 
(fisogastra) of the posterior abdomen is apparently a process of extreme 
difficulty for the virgin female.® In these isolations the females which 
hatched in the culture were very weak, and swelling was achieved by only 
two of them. These two produced a large number of eggs, and all of 
their 612 determinable offspring were males. Contrary to Reuter’s nega- 
tive findings,® it is concluded from these data that Pediculopsis graminum 
(Reut.) is facultatively parthenogenetic. Virgin females give rise only 
to males. 

Reuter, basing his observations on studies of the egg cleavages as well 
as later somatic divisions, states that the diploid chromosome number is 
four.‘ The primary odcytes in the stocks under investigation, however, 
show three tetrads at diakinesis, while the karyokinesis of the primary 
oécyte is invariably characterized by the division of three tetrads (Fig. 1, a). 
Early cleavages of the fertilized egg and somatic mitoses of the advanced 
embryo present complements of six chromosomes (Fig. 1, d, c). Either 
Reuter erred in his determination of the diploid chromosome number of 
Pediculopsis graminum (Reut.), or, contrary to his opinion, the North 
American form here under investigation actually represents a different 
species with a greater chromosome number. If the latter suggestion 
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actually be the case, then the name Pediculopsis dianthophilus (Wolcott) 
should be applied to the American species. 

It has been pointed out above that males represent approximately 5% 
of the normal brood. A count of the chromosome numbers of 180 embryos 
chosen at random showed 9, or 5%, of these to have haploid sets of three 
chromosomes (Fig. 1, 6). Furthermore, a female containing but a single 
embryo was fixed. This embryo was in an early cleavage stage, and its 
sex was not determinable from its gross appearance. Microscopic in- 
vestigation showed this embryo to possess the haploid complement of 
three chromosomes. It has been demonstrated above for cases in which 
females bear but one offspring that the unique offspring are male. From 
this statistical correspondence of the cytological data with those of 
the sex ratio, it seems highly probable that the males are haploid, and that 
Pediculopsis graminum (Reut.) is arrhenotokous. 

This investigation has been carried on in the laboratory of Professor 
Franz Schrader, to whom I make grateful acknowledgment. 


1 Cheyletus eruditus Schr. and Amblyomma dissimile Koch are quite generally believed 
to be firmly established as parthenogenetic. However, the possibility that these species 
may be self-fertilizing hermaphrodites seems neither to have been ruled out nor con- 
sidered. 

2 Schrader, F., Arch. Mikr. Anat., 97, 610 (1928). 

3 Patau, K., Naturwiss. 22, 648 (1934); Zool. Jahrb. 56, 277 (1936). 

4 Reuter, E., Acta Soc. Sct. Fenn., 37, 3 (1909). 

5 Reuter, E., Acta Soc. Sci. Fenn., 36, 1 (1909). See also, Stewart, F. G., and H. E. 
Hodgkiss, N. Y. Agr. Exp. Sta., Tech. Bull., 7, 83 (1908). 

6 Reuter. E., Acta Soc. Sci. Fenn., 36, 1 (1909). 

7 Bradbury, O. C., Jour. Elisha Mitchell Sci. Soc., 42, 94 (1926). 

8 Brucker, C. A., Bull. Sci. France et Belg., 35, 365 (1901). 

® Monterosso, B., Riv. Biol., 16, 80 (1934). 


THE GENERA OF THE ARTICULATED CORALLINES 
By ARTEMIO VALDERRAMA MANZA 
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Communicated December 21, 1936 


In 1904, Dr. Anna Weber-van Bosse (Siboga Exp., Monogr. 61, pp. 
78-110) made certain revisions of the generic limits of the Articulated 
Corallines which direct the way toward a more consistent codrdination 
of both the vegetative and the reproductive characters from microscopic 
as well as macroscopic viewpoints. In connection with certain work on 
the Articulated Corallines collected from the Pacific Coast of North 
America, it has been found necessary to extend the Weber-van Bosse 
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system and to make even more clear the lines of demarcation between the 
different genera recognized. This has been accomplished by instituting 
certain additional combinations of characters both macroscopic and micro- 
scopic, such as modes of: branching, shapes of the segments, positions of 
the conceptacles, positions of the conceptacular pores and the varying 
microscopic structure in both the genicula and the intergenicula. In 
consequence, the following tentative outline is set forth: 


A. Fronds with conceptacles all lateral, 
I. Conceptacles scattered over the intergenicula, 

a. Genicula unizonal or multizonal; intergenicular medullary 
filaments straight, with cells in long and short transverse 
zones, 

1. Amphiroa Lamx. (emend. Weber-van Bosse). 

Lamouroux, Nouv. Bull. des Sci. par la Soc. Philomat., 3, p. 186, 1812; 
Hist. Polyp. Flex., pp. 294-302, 1816; Weber-van Bosse, Siboga Exp., 
Monogr. 61, pp. 86-101, 1904. 

Type of the genus: Amphiroa fragilissima (L.) Lamouroux, Hist. Polyp. 

Flex., p. 298, 1816. 

Type locality: ‘‘O. Americano.” 
b. Genicula multizonal; intergenicular medullary filaments straight, 
with cells in transverse zones of uniform length, 
2. Metagoniolithon Weber-van Bosse, Siboga Exp., Monogr. 61, pp. 86, 

101-104, 1904. 

Type of the genus: Metagoniolithon charoides (Lamx.) Weber-van Bosse, 

Siboga Exp., Monogr. 61, p. 102, 1904. 

Type locality: ‘‘Australasie.’’ 
c. Genicula unizonal, 
aa. Intergenicular medulla unizonal, 
3. Lithothrix Gray, Jour. Bot. 5, p. 33, 1867. 
Type of the genus: Lithothrix Aspergillum Gray, Jour. Bot. 5, p. 33, 
figs. a, b, 1867. 
Type locality: ‘‘Vancouver’s Island.” 
bb. Intergenicular medulla multizonal, 
1. Medullary filaments straight, with cells in transverse 
zones of uniform length, 
la. Segments cylindrical, 
4. Pachyarthron gen. nov. 

Frondes fragiles, dichotomae aut irregulari-dichotomae; segmentis 
cylindricis aut lente compressis; geniculis unizonalibus; filamentis medul- 
laribus intergeniculorum rectis, cellulis in zonis transversis longe equalibus 
positis; conceptaculis semi-globosis in articulis intergeniculorum sparsis. 
Type of the genus: Pachyarthron cretaceum (P. et R.) comb. nov. Coral- 

lina cretacea P. et R., Ill. 40, p. 20, fig. 104, 1840. 
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Type locality: ‘‘ad Insulam Unalaschka.” 
1b. Segments compressed, 


5. Bossea gen. nov. 

Frondes fragiles, dichotomae, pinnatae aut pinnato-dichotomae; seg- 
mentis in parte inferiore compressis aut cylindricis, in parte superiore 
compressis; geniculis unizonalibus; filamentis medullaribus intergenicu- 
lorum rectis, cellulis in zonis transversis longe equalibus positis; concepta- 
culis semi-globosis usque ad conicis in superficiebus planis articulorum 
intergeniculorum sparsis. 

Type of the genus: Bossea plumosa sp. nov. 

Frondibus erectis, laete purpureis, 3-6 cm. longis, oppositis plumoso- 
pinnatis usque ad bi-tripinnatis, infero ramellis curtissimis, supero ramis 
longis uni-bipinnatis, iisdem infero ramellis curtissimis et supero ramulis, 
vestitis; conceptaculis conicis 2-4 in intergeniculo quoque, vulgo in super- 
ficie quaque 2, singulo aut binis (in seriebus longitudinalibus per axes), 
et singulis et centralibus in superficiebus planis intergeniculorum ramel- 
lorum, poris mediis et tetrasporangiis basalibus. 

Species type: tetrasporic, Moss Beach, Coast of Central California, 

North America (Herb. Univ. Calif. No. 545710). 

2. Medullary filaments flexuous and interlacing, 


6. Calliarthron gen. nov. 

Frondes fragilissimae, dichotomae, pinnatae, aut dichotomo-pinnatae; 
segmentis cylindricis aut compressis; geniculis unizonalibus; filamentis 
medullaribus intergeniculorum flexuosis intertextisque; conceptaculis 
semi-globosis aut conicis in superficiebus intergeniculorum, aut in mar- 
ginibus lateralibus aut etiam in superficiebus planis positis. 

Type of the genus: Calliarthron cheilosporioides sp. nov. 

Frondibus 10-29 cm. altis, opposite et laxe pinnatis aut subalterne 
ramosis; conceptaculis in lobis superis et marginibus lateralibus primo 
1-3 in margine quaque oriendis, ultimo quoque in superficiebus planis 
intergeniculorum congestis, conicis, in plantis tetrasporangiiferis apicibus 
obtusis, in plantis antheridiiferis apicibus acutis. 

Species type: tetrasporic, Pacific Grove, Coast of Central California 
(Herb. Univ. Calif. No. 545724); antheridial, Pebble Beach, Carmel 
Bay, Coast of Central California, North America (Herb. Univ. Calif. 
No. 545721). 

II. Conceptacles on the upper lobes of the intergenicula, 


7. Cheilosporum Areschoug, in J. G. Agardh, Species Algarum, 2, sect. 2: 
pp. 543-547, 1852. 

Type of the genus: Cheilosporum sagitiatum (Lamx.) Areschoug, in J. G. 
Agardh Species Algarum, 2, sect. 2: p. 545, 1852. 

Type locality: ‘‘Ile-de France.” 
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B. Fronds with terminal conceptacles, 

I. Fronds with both terminal and lateral conceptacles, 

8. Joculator gen. nov. 

Frondes fragiles, pinnatae usque ad pinnato-decompositae; segmentis 
infero compressis aut cylindricis, supero compressis; geniculis unizonali- 
bus; filamentis medullaribus intergeniculorum rectis, cellulis in zonis 
transversis equalibus positis; conceptaculis terminalibus in ramulorum 
apicibus et etiam in superficiebus planis intergeniculorum lateralibus. 
Type of the genus: Joculator pinnatifolius sp. nov. 

Frondibus erectis, gracilibus, atropurpureis, 4-7 cm. longis, oppositis 
et plumoso-pinnatis, aliquando tripinnatis; conceptaculis in apicibus 
ramulorum aliquando antenniferis, aliquando non-antenniferis et etiam 
superficiebus intergeniculorum in utraque superficie 2, quaque singula 
prope lobam superiorem positis. 

Species type: tetrasporic, Orange County, Coast of Central California, 
North America (Herb. Univ. Calif. No. 545769). 

II. Fronds with conceptacles all terminal, 
a. Branching strictly dichotomous, 

9. Jania Lamouroux, Nouv. Bull. des Sci. par la Soc. Philomat. 3, p. 186, 
1812; Hist. Polyp. Flex., pp. 266-274, 1816. 

Type of the genus: Jania rubens (L.) Lamouroux, Nouv. Bull. des Sci. 
par la Soc. Philomat., 3, p. 186, 1812; Hist. Polyp. Flex., pp. 271-273, 
1816. 

_ Type locality: “‘O. Europaeo.” 

b. Branching wholly or partially pinnate, 
aa. Conceptacles non-antenniferous, 

10. Corallina L. (emend. Lamx.), Linnaeus, Syst. Nat. ed. 12, 1: part 2, 
pp. 1304-1306, 1767; Lamouroux, Nouv. Bull. des Sci. par la Soc. 
Philomat. 3, pp. 185-186, 1812; Hist. Polyp. Flex., pp. 275-292, 1816. 

Type of the genus: Corallina officinalis Linnaeus, Syst. Nat., ed. 12, 1: 
part 2, p. 1304, 1767. 

Type locality: ‘‘O. Europaeo, Mediterraneo.” 

bb. Conceptacles antenniferous, 

11. Corallina (Cornicularia) subgen. nov. 

Frondes fragiles, in toto aut ex parte pinnatae; segmentis cylindricis 
aut compressis; geniculis unizonalibus; filamentis medullaribus inter- 
geniculorum rectis, in zonis transversis equalibus positis; conceptaculis 
terminalibus, omnino aut plerumque antenniferis. 

Type of the subgenus: Corallina (Cornicularia) gracilis Lamouroux, 
Hist. Polyp. Flex., p. 288, pl. 10, fig. 1 (a and b), 1816. 

Type locality: ‘‘Australasie.”’ 

c. Branching pinnate-cymoid, 
aa. Conceptacular pores apical, 
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12. Arthrocardia Areschoug, in J. G. Agardh Species Algarum. 2, sect. 2: 
pp. 547-553, 1852. 

Type of the genus: Arthrocardia corymbosa (Lamk.) Areschoug, in J. G. 
Agardh Species Algarum, 2, sect. 2: pp. 550-551, 1852. 

Type locality: ‘‘ad littora Americae,”’ fide Lamarck; “‘ad oram Capensem 
praesertim in sinu Tabulari nec non Algoensi non infrequenter,’’ fide 
Areschoug. 

bb. Conceptacular pores slightly lateral (slightly below the 
apices), 

13. Duthiea gen. nov. 

Frondes fragiles, pinnato-cymoideae; segmentis basi compressis aut 
cylindricis, supero compressis; geniculis unizonalibus; filamentis medul- 
laribus intergeniculorum rectis, cellulis in zonis transversis equalibus 
positis; conceptaculis in ordine cymoideo-terminalibus, poris oblique- 
lateralibus. 

Type of the genus: Duthiea Setchellii sp. nov. 

Frondibus erectis, 4-6 cm. altis, parte vegetabili opposite pinnatis aut 
tripinnatis, fructifera cymoideo-ramosis; conceptaculis apicalibus in 
intergeniculis terminalibus, intergeniculis conceptaculiferis e marginibus 
superioribus intergeniculorum fertilium in ordine, more cymarum oriendis, 
poris lente oblique-lateralibus, tetrasporangiis in fundo et lateribus cav- 
ernarum conceptaculorum positis. 

Species type: tetrasporic, Blauwklip, near the mouth of Grote River, 
South Africa (Herb. Univ. Calif. No. 545765; Duthie No. 8012). 


AN EOCENE TITANOTHERE FROM SAN DIEGO COUNTY, 
CALIFORNIA, WITH REMARKS ON THE AGE OF THE POWAY 
CONGLOMERATE 


By CHESTER STOCK 
BALCH GRADUATE SCHOOL GEOLOGICAL SCIENCES, CALIFORNIA INSTITUTE TECHNOLOGY 


Communicated January 16, 1937 


Introduction.—The occurrence of Eocene vertebrates in the Poway 
Conglomerate of San Diego County, California, was first brought to my 
attention several years ago when, in the course of excavations north of 
La Mesa and approximately 10 miles east and slightly north of the city of 
San Diego, fragmentary remains of vertebrates were encountered at a 
depth of 40 feet. Examination of this material, permitted through 
courtesy of Director Clinton G. Abbott of the San Diego Museum of 
Natural History, elucidated the information that several mammals found 
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at this locality are either closely related to or identical with forms recog- 
nized in the upper Eocene fauna from the Sespe deposits, north of the Simi 
Valley, Ventura County. 

Further field investigations conducted in the area have uncovered addi- 
tional fossil vertebrate localities and at one of these, during the past 
summer, was found the titanothere to be described in this paper. The 
specimen is noteworthy because it permits a comparison with related types 
from the North American Eocene, so comprehensively discussed by 
Osborn! in his monograph on the Titanotheres. In addition, the presence of 
this material in the Poway Conglomerate, a deposit containing marine 
lenses from which Hanna? obtained an invertebrate fauna tentatively cor- 
related with the Tejon Eocene, offers excellent opportunity to compare the 
age determination of the Poway based on a study of the titanothere with 
that derived from a study of the marine mollusks. 


DESCRIPTION OF MATERIAL 
SUBFAMILY DOLICHORHINAE 
Metarhinus(?) pater, n. sp. 


Type Specimen.—The right half of the muzzle of a skull from the an- 
terior part of the orbit to the anterior end of the maxillary with the canine 
and P1-M3; No. 2037 Calif. Inst. Tech. Coll. Vert. Pale.; Plate 1. Speci- 
men collected by Jack Dougherty. 

Locality.—Sandstones of the Poway formation, exposed in west bank of 
the San Diego River, approximately one-quarter mile north and east of 
the Mission San Diego; Calif. Inst. Tech. Vert. Pale. Loc. 249. 

Specific Characters—Muzzle elongate; naso-maxillary recess deep and 
reaching back of antorbital foramen; antorbital foramen large and situated 
above anterior border of M2. P1 and P2 small; P4 relatively long in 
comparison to its width. Size larger than Metarhinus fluviatilis and 
M. riparius, similar to that of M. earlet and smaller than Mesatirhinus 
superior. : 

Comparison with Metarhinus.—The principal point of difference between 
No. 2037 and skulls of Metarhinus is presented in the great backward 
extent of the naso-maxillary notch as seen in side view. The ascent of the 
lower border of this notch in its posterior course is more gradual than in 
skulls of Metarhinus with which comparison has been made and, moreover, 
extends farther behind the level of the antorbital foramen than in these 
specimens. Above and to the inner side of the foramen is a broad plat- 
form of bone which is grooved longitudinally. This forms the posterior 
part of the floor of the recessed area, the roof of which was formed by the 
nasal bones, unfortunately not preserved in the Poway specimen. A 
similar platform is seen in the skull of Metarhinus fluviatilis, but the lateral 
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recess extends back to a point above and in front of the antorbital foramen. 
The lateral wall of the maxillary in front of the antorbital foramen shows no 
indentation or groove. The medial dorsal contact of the premaxillaries is 
long in fore and ait direction. 

A deep naso-maxillary recess occurs.in the skull of M. riparius, the recess 
reaching nearly to the level of the inner anterior border of the orbit. 
No. 2037 from the Poway differs principally in greater elongation of 
snout as expressed by greater length in front of orbit and by a longer 





PLATE 1 


Metarhinus (?) pater, n. sp. 


Type specimen, No. 2037, Calif. Inst. Tech. Coll. Vert. Pale. Ventral and lateral 
views, X 1/2. 
Poway Conglomerate, San Diego Co., California. Upper Eocene. 


diastema between P1 and canine. In M. riparius the lateral face of the 
snout curves inward more gradually in its upward course and a small plat- 
form extends forward from the antorbital foramen. 

The surface of the maxillary below the orbit in No. 2037 shows that the 
jugal has been broken away along the plane of the suture. Examination of 
the specimen leads to the impression that the infra-orbital shelf was 
fully as well developed as in Metarhinus and thus differed from that in 
Rhadinorhinus. 
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No. 2037 resembles species of Metarhinus in the close approach made by 
the antorbital foramen to the orbital border, in the position of the post- 
narial notch opposite or slightly in advance of the posterior level of M2, 
and in the small size of Pl. The San Diego specimen differs in presence of 
a deeper naso-maxillary recess, longer snout, longer postcanine diastema, 
greater anteroposterior diameter of P4 in relation to its transverse diameter, 
and in a more posterior position of the antorbital foramen since the poste- 
rior rim of this opening is situated above the anterior end of M2. A well 
developed hypocone is present in M3. 

Comparison with Mesatirhinus.—The tooth-row in No. 2037 resembles 
in length that in the type of Mesatirhinus petersoni and in the type of M. 
superior, but the premolar series is shorter. The canine is larger than the 
comparable tooth in M. superior but may be similar in size to that in 
M. petersoni. P1 is smaller in the San Diego specimen than in Mesa- 
tirhinus, but P2 has the subtriangular shape seen in this genus. 

In P3 and in P4 the internal cingulum swings around from the anterior 
and posterior sides of the tooth, but does not extend entirely along the 
inner side of the base of the principal internal cusp. M2 in No. 2037 is 
narrower relative to its length than is the comparable tooth in Mesatirhinus. 
The postnarial notch reaches slightly farther forward than in M. superior, 
but not quite so far as in M. petersoni. The antorbital foramen is large, 
relatively more so than in Mesatirhinus. The distance between the lateral 
border of this opening and the anterior border of the orbit is shorter than in 
Mesatirhinus. The narial opening is much more deeply recessed as seen 
from the side than in the latter genus. 

Comparison with Dolichorhinus——No. 2037 is distinctly smaller than 
typical members of the genus Dolichorhinus and differs markedly from these 
in the much more deeply recessed narial opening, absence of a secondary 
palate, more posterior position of the antorbital foramen, decidedly shorter 
distance between this foramen and the orbital border, shorter postcanine 
diastema and in the subtriangular rather than subquadrate shape of P2. 
In specimens of Dolichorhinus the cingulum is complete along the inner 
sides of the posterior premolars, not interrupted as in No. 2037. 

Comparison with Rhadinorhinus.—The tooth-row in No. 2037 resembles 
in length that in the type of R. diploconus, but lacks the upward curvature 
in a vertical plane of the anterior cheek-teeth. P2 is not subquadrate as 
in Rhadinorhinus; the postcanine diastema is longer; the postnarial notch 
is farther back; the infra-orbital shelf is better developed in the San 
Diego specimen. While the lateral recess of the nose is deep in Rhadi- 
norhinus, the notch does not extend backward so far as in No. 2037 and 
the lower posterior rim of the notch is sharper. The antorbital foramen is 
large as in Rhadinorhinus and, as in this genus, is closely situated to the 
border of the orbit, but has a position farther back with reference to M1. 
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MEASUREMENTS (IN MILLIMETERS) OF No. 2037 


Length from anterior end of C to posterior end of M3 200 
Length of premolar-molar series through middle 166.4 
Length of molar series through middle 100.7 
Length of postcanine diastema 14.4 
C, anteroposterior diameter 20 
C, transverse diameter 17.5 
P1, anteroposterior diameter 12.4 
P1, transverse diameter 7.2 
P2, anteroposterior diameter 18.6 
P2, transverse diameter 16.4 
P3, anteroposterior diameter through middle 17.4 
P3, transverse diameter 21 
P4, anteroposterior diameter through middle 20.7 
P4, transverse diameter 25.4 
M1, anteroposterior diameter through middle 27 
M1, transverse diameter 30 
M2, anteroposterior diameter through middle 33 
M2, transverse diameter 32.8 
M3, anteroposterior diameter through middle (approximate) 39 
MB, transverse diameter 39 


Summary of Comparisons.—Among North American Eocene titanotheres, 
the type, No. 2037, resembles most closely the genera Metarhinus, Mesa- 
tirhinus, Dolichorhinus and Rhadinorhinus. In its most striking character 
the Poway specimen resembles more closely the first and last genera than 
it does Mesatirhinus. Resemblance to Dolichorhinus is even more remote. 
No. 2037 represents a skull with elongated muzzle and with very deeply re- 
cessed naso-maxillary opening. In latter character the San Diego species 
is distinctly more modified than any known member of the Dolichorhinae. 
Nearest approach in the character of deep recess is to Metarhinus. In 
Rhadinorhinus the opening is also deeply recessed, but in other characters 
of skull and dentition this genus resembles the San Diego form less closely 
than does Metarhinus. The characters which distinguish No. 2037 from 
Metarhinus are possibly diagnostic of a new genus, but I hesitate to de- 
scribe the Poway specimen as such in view of the present fragmentary 
material. No. 2037 may represent an offshoot of the metarhine group, 
later and more specialized than Metarhinus, in which an advance is shown 
particularly by the elongation and further modification of the snout and 
likewise perhaps by the well developed hypocone in M3. 


REMARKS ON THE AGE OF THE PowAy CONGLOMERATE 


Hanna?’ in 1926 published the results of his geologic study of the La 
Jolla Quadrangle. Rocks of Eocene age exposed in this area were recog- 
nized by Hanna as comprising, in ascending order, the following: (1) 
Delmar Sand, (2) Torrey Sand, (3) Rose Canyon Shale and (4) Poway 
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Conglomerate. The three lowermost beds were grouped together as the 
La Jolla formation and were considered as intermediate in age between 
the Meganos and Tejon formations of the Pacific Coast marine Eocene. 
The La Jolla formation was regarded as an apparent correlative of the 
Domengine Eocene of the Simi Valley section, California. 

The locality of occurrence of Metarhinus(?) pater is in an area mapped by 
Hanna as the Poway Conglomerate. The specimen was found in a sand- 
stone stratum between conglomerate layers. In reaching an age deter- 
mination of this portion of the Poway it is worthy of note that the titano- 
theres with which the Poway species is compared, include types that range 
upward stratigraphically from the Metarhinus beds of the Bridger upper 
middle Eocene to and into the Dolichorhinus beds representing the Uinta B 
stage of the upper Eocene, as these horizons are now recognized in the 
Tertiary basins of the Cordilleran Province. Relationship of the San 
Diego species appears to be with the metarhine titanotheres, but the 
former is more specialized and advanced than typical Metarhinus. M.(?) 
pater differs widely from the genera Diplacodon, Eotitanotherium and 
Protitanotherium recorded from the horizon Uinta C of the Uinta basin. It 
is distinctly less advanced and clearly an earlier type of titanothere than 
Teleodus californicust from the uppermost Eocene Sespe of Ventura 
County, California. Thus the relationships of the San Diego titanothere 
suggest an age determination of the land-laid Poway as not older than the 
Bridger upper middle Eocene and not as late as the Uinta C stage of the 
Upper Eocene. On the other hand, the advance in structural characters 
which M.(?) pater exhibits in a comparison to the typical Metarhinus 
may indicate a stage of development later in time than that represented 
by the related forms from Uinta B. Thus it is possible that M.(?) pater 
occurred at a time when Diplacodon and Protitanotherium were present 
during the accumulation of the later Uinta. In either case, the San Diego 
titanothere points to an upper rather than middle Eocene age for the land- 
laid Poway. 


1 Osborn, H. F., U. S. Geol. Surv. Mon., 55, 1, 383-434 (1930). 

2 Hanna, M. A., Univ. Calif. Publ., Bull. Dept. Geol. Sci., 16, 256-263 (1927). 

3 Hanna, M. A., Univ. Calif. Publ., Bull. Dept. Geol. Sci., 16, 187-246, pls. 17-23, 
map (1926). 

4 Stock, C., Proc. Nat. Acad. Sct., 21, 456-462 (1935). 








54 GENETICS: A. G. STEINBERG Proc. N. A. S. 


RELATIONS BETWEEN CHROMOSOME SIZE AND EFFECTS OF 
INVERSIONS ON CROSSING-OVER IN DROSOPHILA 
MELANOGASTER 


By ArtTHUR G. STEINBERG 
DEPARTMENT OF GENETICS, COLUMBIA UNIVERSITY 


Communicated January 15, 1937 


Schultz and Redfield’? have shown that inversions in chromosomes I 
and II cause an increase in crossing-over in chromosome III and that in- 
versions in chromosomes I and III have similar effect on crossing-over in 
chromosome II. Bridges* showed that inversions in chromosomes I, II 
and III when present simultaneously have no effect on crossing-over in 
chromosome IV. Steinberg® showed that inversions in chromosomes II 
and III cause an increase in crossing-over in chromosome I. It was 
found both by Schultz and Redfield and by Steinberg that different in- 
versions caused different amounts of increase in crossing-over in a given 
chromosome. 

Table 1, which is a combination of data of Schultz and Redfield (un- 
published) for chromosomes II and III and Steinberg’ for chromosome I, 
shows that a given inversion causes a different average increment in 
crossing-over per unit map length in different chromosomes. For ex- 
ample, the C1B inversion causes a thirty-eight per cent increase in crossing- 
over per unit map length in chromosome II while in chromosome III it 
causes an increase of forty-six per cent per unit map length (table 1). The 
question arises as to what causes this differential quantitative effect of 
inversions on crossing-over in non-homologous chromosomes. While no 
complete answer to this question can be given at the present time, refer- 
ence to tables 2 and 3 will show a partial answer. 


TABLE 1 


SHOWING A SUMMARY OF THE INCREASES IN CROSSING-OVER CAUSED IN Non-HoMOoLo- 
GOoUS CHROMOSOMES BY INVERSIONS (VALUES EXPRESSED AS DECIMAL FRACTIONS) 


CHROMOSOME CHROMOSOME CHROMOSOME 
I II III 


INVERSIONS TOTALS 
C1B ply 0.38 0.46 0.84 
Curly 0.22 cat 0.49 0.71 
Payne 0.33 0.55 peer 0.88 


Table 2 shows a comparison of the ratios of the relative lengths of meta- 
phase chromosomes‘ with the ratios of the interchromosomal effects 
of an inversion or pair of inversions on crossing-over in these chromosomes. 
Thus chromosome I is sixty-nine per cent as long as chromosome II and 
the effect of the Payne inversions on I is sixty per cent of the effect of the 
same inversions on II, likewise chromosome I is fifty-seven per cent as 
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long as III while the effect of the Curly inversions on the former is forty- 
five per cent of that on the latter; and finally chromosome II is eighty- 
three per cent of III and the effect of the C1B inversion on chromosome II 
is also eighty-three per cent of that on III. Within the limits of error which 
are probably large because: 

(1) of the difficulties involved in measuring Drosophila chromosomes, 

(2) the interchromosomal effects of the inversions were measured in 
three different experiments, and 

(3) the first chromosome was not completely marked, 
these pairs of ratios may be considered very similar. 


TABLE 2 


SHOWING THE RATIOS OF THE RELATIVE METAPHASE LENGTHS OF THE GONIAL CHROMO- 
SOMES AND THE CORRESPONDING RATIOS OF THE INTERCHROMOSOMAL EFFECT OF 


INVERSIONS 
CHROMOSOMES . We L/III II/Ill 
Ratio of sizes 100/145 = 0.69 100/175 = 0.57 145/175 = 0.83 
Ratio of effects 0.33/0.55 = 0.60 0.22/0.49 = 0.45 0.38/0.46 = 0.83 


TABLE 3 


SHOWING THE PERCENTAGE OF THE TOTAL CHROMATIN LENGTH CONTAINED IN EACH 
OF THE CHROMOSOMES AND THE CORRESPONDING PERCENTAGE OF THE TOTAL INTER- 
CHROMOSOMAL EFFECT OF INVERSIONS SHOWN BY EACH CHROMOSOME 


CHROMOSOMES I II III IV 
Percentage of total chromatin 23.3 33.7 40.7 2.3 
Percentage of total effect 22.6 38.3 39.1 0.0 


Table 3 shows a comparison of the percentage length of a given chromo- 
some (i.e., the length of the chromosome divided by the sum of the lengths 
of all four chromosomes and multiplied by 100) with the per cent of the 
total increase in crossing-over shown in this chromosome. For example, 
the increases in crossing-over in chromosome I are .22 + .33 = .55 (table 1), 
the total increase in crossing-over is .84 + .71 + .88 = 2.43 (table 1); 
the percentage of the total effect which is shown by chromosome I is 
.55/2.43 X 100 = 22.6. These figures (table 3) show even more clearly 
than those in table 2 that the effect of an inversion on crossing-over in a 
non-homologous chromosome is directly proportional to the length of the 
chromosome affected. In addition table 3 shows that the proportion of the 
total effect which goes to a given chromosome is approximately equal to 
that proportion of the total length of all the chromosomes which the affected 
chromosome constitutes. Therefore, the effect of an inversion on crossing- 
over in a non-homologous chromosome is a function of both the in- 
version itself and of the length of the chromosome affected. This fact ex- 
plains Bridges’* observation that inversions do not affect crossing-over in 
chromosome IV. Chromosome IV, according to the above, would be ex- 
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pected to show about two per cent of the total increase in crossing-over 
(table 3), i.e., crossing-over in chromosome IV should be increased by ap- 
proximately five per cent. Such an increase in crossing-over in chromo- 
some IV is so small as to be detectable only in extremely large test popula- 
tions. 

The author wishes to acknowledge his indebtedness to Professors L. C. 
Dunn and D. E. Lancefield for their helpful suggestions, and to Dr. J. 
Schultz and Miss H. Redfield for kindly permitting him to use their un- 
published data. 


1 Morgan, T. H., Bridges, C. B., and Schultz, J., Yearb. Carnegie Instn., 31, 303 


(1932). 

2? Morgan, T. H., Bridges, C. B., and Schultz, J., Yearb. Carnegie Instn., 32, 298 
(1933). 

3’ Morgan, T. H., Bridges, C. B., and Schultz, J., Yearb. Carnegie Instn., 33, 274 
(1934). 


4 Morgan, T. H., Bridges, C. B., and Sturtevant, A. H., Bibliog. Genet., 1 (1925). 
5 Steinberg, A. G., Genetics, 21, 615 (1936). 


LINKAGE STUDIES OF THE RAT (RATTUS NORVEGICUS). II‘ 
By HELEN DEAN KING AND W. E. CASTLE 


WISTAR INSTITUTE AND UNIVERSITY OF CALIFORNIA 


Communicated January 13, 1937 


King? has recently described the occurrence and inheritance of waltzing 
in an inbred race of rats reared at the Wistar Institute. Concerning the 
mode of inheritance this statement was made: ‘Available evidence seems 
to indicate that waltzing is due to a single gene, but that the extent to 
which this gene expresses itself depends upon modifying factors which 
favor or inhibit its action.’’ Selection of the most pronounced waltzers 
through a series of generations increased the incidence of waltzing. In the 
F, generation the percentage of waltzers produced by waltzing parents 
had risen to 93 per cent. 

The waltzing mutation made its first appearance in albino individuals 
and though it has since been introduced into colored individuals, it con- 
sistently occurs in a higher percentage of the albino than of the colored 
individuals in litters containing both sorts. This indicates that it is 
linked with albinism and that its gene is borne in the same chromosome 
as the albino gene. The evidence for this conclusion is as follows: 

Crosses were made between albino waltzers and colored non-waltzers of 
the varieties known as hairless, kinky and blue. Ali F, individuals were, 
as expected, colored non-waltzers. The F: population obtained from each 
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cross is classified in table 1 as regards the occurrence of albinism. No indi- 
cation was found of linkage between waltzing and the other mutant char- 
acters involved in these crosses, hairless, kinky and blue. 


TABLE 1 


F, POPULATIONS IN WuicH Is SHOWN THE COMPARATIVE INCIDENCE OF WALTZING 
AMONG ALBINO AND NON-ALBINO INDIVIDUALS, RESPECTIVELY 
ALBINO COLORED 


NON- ALBINO PER CENT NON- COLORED PER CENT 
WALTZERS WALTZERS WALTZERS WALTZERS WALTZERS WALTZERS 


Cross 1. Albino waltzer X 


hairless 82 7 7.8 202 3 1.5 
Cross 2. Albino waltzer X 

kinky 36 8 19.0 154 10 6.1 
Cross 3. Albino waltzer < 

blue 88 13 14.7 316 22 6.5 


It will be observed that the incidence of waltzing was in all three of the 
F, populations higher among albinos than among their colored sibs. 

The same relation was, with a single exception, found to hold in cer- 
tain back-cross populations produced by mating F; colored non-waltzers 
with albino waltzers, as shown in table 2. 


TABLE 2 


COMPARATIVE INCIDENCE OF WALTZING IN ALBINO AND COLORED INDIVIDUALS IN 
Back-Cross POPULATIONS 


ALBINO COLORED 
NON- ALBINO PER CENT NON- COLORED PER CENT 
WALTZERS WALTZERS WALTZERS WALTZERS WALTZERS WALTZER3S 


Cross 4A. Albino waltzer < 


Curly 104 21 16.8 117 21 15.2 
Cross 4B. Albino waltzer X 

Curly 24 18 42.8 35 9 20.4 
Cross 5A. Albino waltzer X 

Black 44 16 26.6 64 16 20.0 
Cross 5B. Albino waltzer X 

Black 58 10 14.7 63 12 16.0 
Cross 7. Albino waltzer X 

Curly: 52 38 42.2 47 24 33.8 

Totals 282 103 26.7 326 82 20.1 


In the Cross 5B population alone was the percentage of waltzers as high 
among the colored as among the albino individuals. For all the popula- 
tions of table 2 combined, 26.7 per cent of the albinos were waltzers, but 
of the colored individuals only 20.1 per cent. 

One further point of interest should be noted in connection with Cross 
3, table 1. The blue parent to this cross was heterozygous for yellow, a 
recessive gene borne in the same chromosome as the albino gene. If 
waltzing is linked with albinism, it should show linkage also with yellow. 
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The F: population derived from Cross 3 may be separated into three groups 
as shown in table 3. 
TABLE 3 


A FuRTHER CLASSIFICATION OF THE F; POPULATION OF 426 INDIVIDUALS DERIVED FROM 
Cross 3, TABLE 1 


NON- PER CENT 
WALTZERS WALTZERS WALTZERS 
Gray, black or blue 
(i.e., non-yellow colored) 259 16 5.8 
Yellow or cream (dilute yellow) 57 6 9.5 
Albino | 75 13 14.7 
Totals 391 35 8.2 


As expected, the percentage of waltzers in the yellow group is higher than 
that in the non-yellow group, suggesting linkage of waltzing with yellow 
as well as with albinism. 

The data presented in tables 1-3 had been obtained prior to March 17, 
1936. Subsequent experiments support the tentative conclusion reached 
at that time, that linkage exists between waltzing and albinism. See 
table 4. : 

TABLE 4 
FURTHER DATA ON THE COMPARATIVE INCIDENCE OF WALTZING AMONG COLORED AND 
ALBINO Sips IN Back-Cross POPULATIONS 
ALBINO COLORED 


NON- ALBINO PER CENT NON- COLORED PER CENT 
WALTZERS WALTZERS WALTZERS WALTZERS WALTZERS WALTZERS 


Cross 7A. Albino waltzer X 


Curly: 171 116 40.4 176 92 34.3 
Cross 5C. F, (Albino waltzer 

X black) X black waltzer 30 9 23.1 167 43 20.5 
Cross 5D. Albino waltzer X : 

black 146 38 20.6 159 38 19.3 
Cross 2A. Albino waltzer X 

kinky 56 21 27.3 68 16 19.0 

Totals 403 184 31.3 570 189 24.9 


In each of the four back-cross populations of table 4, the percentage of 
waltzers is higher among albinos than among their colored sibs, the totals 
being 31.3 per cent among the albinos to 24.9 per cent among the colored. 

Further experiments to test for linkage between waltzing and yellow are 
as yet incomplete and inconclusive. The total back-cross population ob- 
tained by summarizing the various experiments described in tables 2 and 4 
consists of 2139 individuals, which fall into the following classes: 


ALBINO COLORED 
NON- ALBINO NON- COLORED e 
WALTZERS WALTZERS WALTZERS WALTZERS } 
i 
Observed numbers 685 287 896 271 ; 


Corrected for overlapping 463 550 606 519 ‘ 
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When we attempt to calculate the crossover percentage between al- 
binism and waltzing, we are confronted with this difficulty, that the ob- 
served number of waltzers (558) is less than the observed number of non- 
waltzers (1581), with which it should be theoretically equal. This, as al- 
ready explained, we may assume to be due to the occurrence of normal 
overlaps among the non-waltzers, i.e., of individuals genetically waltzers 
but actually not manifesting the waltzing behavior. Nevertheless, we 
observe that the non-crossover class, colored non-waltzers, is larger than 
the crossover class, albino non-waltzers, and also that the other non- 
crossover class, albino waltzers, is larger than the other crossover class, 
colored waltzers, these being the relations which we should expect if link- 
age exits between albinism and waltzing. 

If we assume that overlapping to normality occurs equally among the 
albino and the colored waltzers, we may use the observed frequencies in 
calculating a crossover percentage. 


The observed number of crossovers will accordingly be the sum of the 
albino non-waltzers (685) plus the colored waltzers (271), or 956, which is 
44.7 + 0.7 of the total population, 2139. This is a very loose indicated 
linkage but one of undoubted significance in view of the small probable 
error and the consistency of the occurrence of a larger percentage of 
waltzers among albino than among colored individuals in eight out of nine 
of the back-cross populations. 


This calculation of the crossover percentage may be checked by actually 
deducting from the classes of non-waltzers, and adding to the classes of 
waltzers, enough individuals (considered as normal overlaps) to make the 
total in each of the two groups equal to the half population, 1019. Sucha 
“correction for overlapping’’ has been made in the foregoing summary of the 
backcross populations. If we base a calculation of the crossover percentage 
on these corrected numbers, we get 463 + 519 = 982 as the total of the 
crossover classes, which is 45.8 + 0.7 per cent of the entire population, 
again a statistically significant indication of linkage. 

It has been known for some time that the albino chromosome of the rat 
contains genetic loci for two other mutant characters, red-eyed yellow (r) 
and pink-eyed yellow (p). We are now warranted in assuming the exis- 
tence in this same chromosome of a fourth gene which in its mutant form 
may produce waltzing. The order of the four genes and their approxi- 
mate map-distances are as follows: 


Cc R e W 
0 0.3 20 45 





The experiments described in this paper have been made by King at the 
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Wistar Institute in consultation with Castle, who has formulated this 
report. 


1 See also Proc. Nat. Acad. Sct., 21, 390-399 (1935). 
2 King, H. D., Jour. Mammology, 17, 157-163 (1936). 
3 Castle, W. E., Genetics and Eugenics, 4th Ed., 234 (1930). 


THE COEFFICIENT OF PREDICTION-ACCURACY 


Computation of the Portion of Individual Cases Common to Both the Particular 
Prediction-Distribution and to the Subsequently Determined Actual-Distribu- 
tion of the Same Measured and Counted Phenomena 


By Harry H. LAUGHLIN 


EuGENICS RECORD OFFICE, CARNEGIE INSTITUTION OF WASHINGTON, COLD SPRING 
HARBOR, LONG ISLAND, N. Y. 


Communicated January 4, 1937 


Consider first the specifications for a sound formula for the coefficient of 
prediction-accuracy. Regardless of how, where or why we got it, if we 
have in hand a prediction in terms of a probability-distribution, and if the 
subsequently determined actual-distribution corresponds exactly with the 
prediction-distribution, it is evident that the particular prediction is per- 
fect, and that consequently the coefficient of prediction-accuracy should 
equal 1.000. If, in another case, there be no overlapping of the two proba- 
bility-distributions (that is, if the subsequently determined actual-distri- 
bution does not overlap to any degree with the given prediction-distribution 
which is being tested), then the coefficient of prediction-accuracy should 
work out in the value 0.000. If, in still another case, 95 per cent of the 
individual values actually found were as predicted, then the coefficient of 
prediction-accuracy should register 0.9500. And so on, if three-fourths or 
if one-half of the case values are common between the prediction-distribu- 
tion and the actual-distribution, then the nature of the coefficient of pre- 
diction-accuracy should be such that its formula would produce measured 
coefficients respectively of 0.7500 and of 0.5000. Thus the CPA is prop- 
erly measured by that portion of cases common to both areas “‘p’’ and ‘‘a.”’ 

Also both the particular prediction-basis and the subsequent actual 
check-up by first-hand survey, must be plotted in the same coérdinate 
frame, in which K = f(R), in which R equals the measure of subject- 
quality and K equals the probability-value. Also each distribution-area 
(that is, ‘‘p’’ the prediction-area and ‘‘a’’ the actual area) must equal 1.000. 

Whenever any two normal probability-curves, each of which bounds an 
area which equals 1.000, are plotted within the same coérdinate frame (as 
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are the cases herewith) there will always be some overlapping. If the 
prediction-area be an ‘‘absolute guess or a chance bet”’ such set-up should 
cause the CPA to ‘move away from 1.000 or perfection toward but, of 
course, never reaching zero”’ in value. 


aK 
































K K=f(R) Kfe tKfe 
A Normal Gaussian 
——___ > R fe 
FIGURE 1 
Kye being constant, the area varies directly The Standard Deviation ¢ 
with o, the Standard Deviation. being constant, the area 


varies directly with the Kye. 


When two distribution-areas are plotted within the same frame of co- 
ordinates they produce four points on the R-or-subject-scale which are of 
basic importance in the determination of the coefficient of prediction-accu- 
racy. These four points indicate respectively the positions of 


(a) fc or the fluctuation-center of R-values in ‘‘p’”’ the prediction-distri- 
bution, 

(b) fc or the fluctuation-center of R-values in ‘“‘a’”’ the actual-distribu- 
tion, 

(c) the base of the left-most common probability-value caused by the 
intersection of the boundary curves of the two areas ‘‘p’’ and “a,” and 

(d) the base of the right-most common probability-value caused by the 
intersection of the boundary curves of the same two areas “‘p’”’ and “‘a.”’ 


Given the foregoing specifications of a coefficient of prediction-accuracy, 
we find that the possible situations fall into three cases depending upon the 
relative positions of the four above-listed points on the same R-scale, and 
upon the relative K,. values in “‘p” and ‘‘a’’—that is, in the prediction and 
actual areas. 

With the two ‘‘fc’s’”’ and their respective area-formulae given, it is evi- 
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THE COEFFICIENT OF PREDICTION -ACCURACY (CPA) 


THE PORTION OF CASES COMMON TO BOTH THE PARTICULAR PREDICTION -DISTRIBUTION AND ITS 
SUBSEQUENT ACTUAL-DISTRIBUTION OF MEASURED AND COUNTED PHENOMENA. 


CASE I. 
with Example. 


Each is plotted 


PROPERTIES ae total Prediction-area and the total Actual-area eack +1000. 
within the same codrdinate frame, -the same K, that ts the same ay or Pro- 


bability-scale, and the same R-elass measure of the common subject. 

1 R-class-range = 1 R-unit. 1 Area-unit,= .01 K by 1 R. 
(b) Type of each area-bounding curve, Ke €°?% Fack curve with constant Ke, and 
with tts two O% either equal to or different from each other, but always 


(G+ O)+2 =O symmetrical. 
(¢) Jwo Common Probability -values, due to qrortapping Of curves  # and “a”, one 

cK fo the left of both fcp and fca, the other cK to the right of both fCp and fCa. 
(A) Under the particular probability -curve, whether symmetrical or skewed, 

Area to left fc : Area to right fc =Q:Q when Jotal Area =1.0000. 

Area to the left of ft» —%_- Area to the right of fc= =? 
(e)C.PA*== Four constituent Integral-blocks. Jwo for Area with greater Kye, one in each 
Of its lower extremities, © 7. e and @ Jd °% and two for Area with lesser Kc, in its 
foo right, 


co ye 
central region, @ fA ' ad @ 
Ky hk, 


GIVE iction-Area = 1.000 GIVEN Actual Area = 1.0000 
at point 18.75 on R-scale at point 20.00 onR-scale 
= 0.11250n K-scale = 0.0750 onK-scale 

3.3245 R- 4.9868 R-units 

3.6542 R-units 6.6490 R- units 

R- units 3.3245 R- units 





fe 















fe * fluctuation -center = standard devietion symmetrical  ¢Kz+the left-most probability-value 
common to both p and a. 


LEGEND: K-~ Probability 
R» Measured Subject -quality € = left of particular fe ~ie when probability-curve 1s 
P* @ property of Prediction-eres 4 © right of perticuler fc ota wi sewee cK,» the corresponding right-most 
8+ @ property of Actual-area Kf * probability velue et particular fe Oe + D, = 2.0; when both partic common probability 
sdar area and Kic are constant g *(fep ~fea) in R-units 


FIND 1 = fe predicted ~ fc actual «g = 1250 
(a) R- locus of Ke When (fep ~<Kege) +Xz, ond (fea -CKyg)=Xz,, solve for the positive value of Xz by means of the two simultaneous equations 


Kfep -€ #Kfca -€ een The two formulae for the Common Probability at <Kz, and 
Xz =X_-g Locus Ky on R-scale = locus fe predicted-X,. Xz =4.59 R-units, Xz, = 5.84 R-units 
*. The locus <Ky = 18.750 ~4599 = 14.151 on the R-scale 
(&) R-locus of Ky When (fca ~cKuyn) *X_,and (fep ~cKugh/)*X,,, Solve for the positive value of X, by means of the two Simultaneous equations 
Kfcp -€ = Kfea € HH The two formulae for the Common Probability at cK, , and 
X,,=X,-9 Locus cK, on R-scale = locus fe actual + X,, X_* 1.98 R-units, X,,= 3.23 R-units 
.”, The locus <K, = 20900 + 1.984 = 21.984 on the R-scale 
3 = When the total area under the particular probability-curve = 1.0000 
(a) Predictign-areag t eff of fcp = 2 0, ~ (0p +O) = 1.0992 2 ragos}* #0000 
(b) Prediction-area to right of fcp = 20, ~(Oe+O,)= .QOOB + 2* 4504) ~ 
(c) Actual ~— grea to left of fea = 207 + (0g +0,) = 1.3333 + fr | * 1.0000 
(@) Actual — area to nght of fea = 20,—Op +0) = .6667 + 2*.3334]° ™ 
4. The Value of Constituent Integral-biocks © and @ from the overlapping distribution-orea with the greater Kfc, and of 
Constituent Integral-biocks Q and @ from the overlapping Wistribution-area with tha lesser Kfc 
I a Lf 
(a) When Production area = 10000 Qfcp~cKy+460 R-units » £2985 “Oy —+[5496- (3059 10992] « 1144+ Constituent Intagral-biock @ 
@fcp~cK, +323 R-units = 1.0709 Ay—+[4504-(3599 0.9008) = 1262» Constituent Integrat-biock @ 


(b) When Actual - area = 10000 Ofca~cKy+ 585 R-units> 8796 Ofg———[(3105 1.3393)] > 4140+ Constituent Integrat-biock @ 
@fca~cK,= 1.98 R-umts= $968-0,q———+ [2247 0.6660] _1498*Constituent Integral-bick @ 


£4 Constituent Integral-blocks ———» ,8044» COEFFICIENT of PREDICTION ACCURACY (CPA) 


* 1 *Given equivalent number of Standard Deviations I =Numerical equivalent of Area to the left or right, as the case may be, of the given fc. 


Il-Portion of total symmetrical area subtending given number of Standord Deviations (from Table). 
WV-Total Area,if symmetrical,on basis of particular given J. 
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dent that the first main step toward solution consists in locating the posi- 
tion, on the R-scale, of the foot of each common probability-value. 

The analytic geometry of the problem is developed by the three accom- 
panying plates, one for each of the three cases. In each case the final solu- 
tion is reached by determining the area covered by the overlapping of the 
prediction-probability-distribution-area and the actual-probability-distri- 
bution-area. 

As above stated, the two areas being “‘compared-for-overlapping”’ must 
be plotted in the same frame of codrdinates, and in each case the total 
area under each of the two area-bounding curves must be equal to unity. 
Regardless of the number of cases upon which the particular distribution 
is based, said distribution must be reduced to percentages by R-classes so 
that the total area under the particular bounding-curve is always equal to 


























INOS See Se @ OR Sree eee Kfe K ‘s f(R) 

Normal Gaussian 

O, 
r.6068 ---—-—---------— <— Point of maximum 
K 
>R \ i 
R at one O from fe 
FIGURE 2 


The development of 0.6065 Kye as the K-point at which the Standard Deviation ¢ 
always crosses the K-scale. 


1.000. This rule holds true for both the prediction-distribution and the 
actual-distribution—each must always equal 1.000. Beyond this specifica- 
tion for unity in extent the two areas and their bounding-curves may be 
very different. Each curve may have its own fluctuation-center (fc), its 
own highest probability-value (K,,) and it may be skewed to any degree 
to either the right or the left. 

The three cases, into which the possible types of the present problem 
fall, run as follows: 

Case I is characterized by the set-up in which there are two common 
probability-values. On the R-scale one common probability-value (cK/) 
lies below the two fluctuation-centers (fcp and fca), and the second common 
probability-value (cKr) lies above these two fluctuation-centers. 








64 GENETICS: H. H. LAUGHLIN Proc. N. A. S. 


CASE I. 
with Example. 


OPE (8) Tre t afal Prediction-area and the total Actual-area each =1.0000. _Each is 
d within the same codrdinate frame, ~- the same K, that ts the same crepancy 
or ekadalivaret and the same R-class measure od the common subj 
1 R-class-range = 1 R-unit 1 Area-unit =.01 K by iR. 
(b) Type of each area-bounding curve, K= Kfe £4" Zack curve with constant Kr, and 
with its two O% either equal to or different from each other, but always 
(G+) +2 = 0 symmetrical. 
) Iwo Common Probabiltty-values,due to overlapping of curves “p” and "a", either one 
cK, to the left and the other cK, between *s tt ‘4 C's p and @,or ck between the 
two {c's, and the other cK, to the right of botk 
(a) Under the gaaeee Probability- — whether iad or shewed, 
Area to fC: Area to right {C<G:G when Total Area = 1.0000, 
Area to/left §C =-abq- and Area trope Of fC* abr 
() CPA=£ Tiree Constituent aiepal blocks. Iwo for Area with greater Kfc, in its 
lower extremities, © - as ad @- p lek and one for Area with leser Kfc @ = : iy 


GIVEN Prediction-Ares =1.0000 Gl VEN Actual-Area - 1.0000 
at point 15.00 on R-scale at point 20.00 on R-scale 
= 0.1125 on K- = 0.0750 on K-scale 















= 3.3245 R-units = 4.9868 R-units 
= 2.7448 R-units = 6.6490 R-units 
2 = 3.3245 R-units 














1S) 
°? 
N 
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Sr ed 
eon 
LEGEND: K- Probability fe « fluctuation -center J,= standard deviation symmetrical cKy =the left-most probaotlity - value 
Re Measured Subject-quality J = left of particular fe ~1e. when probability curve is common to beth p and a. 
P ° @ property of Prediction area 4 = right of particular fc te mol beawea” cK, « the corresponding right- most 
; * @ property of Actual-area  Kfce probability value at particulor fc. Og + O, = 205, when both partic common probability. 
vler ora ond Kh. are constont. g » Kp ~fca) in R- units. 


FIND 1. fe predicted ~ fe actual +g = 5.00 units on the R-scale 


2 (a) R-locus of cKy when ep ~cKigt) =X, Gea~Kyg) Xz, , solve for the posttive value of Xz by means of the Simultaneous equations 


gs 
Keep -€ Thee «= Kfea ite The two formulae for the Common Probability af <Kz, and 
Me=Xy,-g. Locus cKy on R-scale = locus fe predicted -Xz. X_=4.71 R-units, Xz, = 9.71 R-untts. 
The locus cKyz = 15.00-4.71 = 10.29 on the R-scale 
w R- locus of cK, when ep ~ cKuge = Xa, and (fca~cKupht)=X_ , solve for the positive value of X, by means of the simultaneous equations: 
Kfep- MOD» Kfea- € . The two formulas for the Common Probability at cK, and 
Xerg-Xze,- Locus cK, on R-scale «locus fe actual ~X,.. X_= 1.39 R-units, X,, = 3.61 R-units. 
-s The locus cK, * 20,00 ~ 1.39 » 18.61 on the R- scale. 
3. When the total area under the particular probability-curre = 1.0000 
(a) Prediction-area to left of fep = 20¢+ (d+ G)« .82563 +2 = 4128)» oo50 
(b) Prediction-area to right of fcp= 20, + (Oz + 0,)*1.17437+2= ‘cont ? 
(c) Actual — area to left of fea 2Q~(Ce+Q) #1.33333+2 « .6667) | dog 
(ad) Actual — grea to right of fca=20,+\O+Q) = .66667+2« peed 
4. The Value, ie. the Area, of Constituent Integrat-biocks ©) ond @ from the overlapping distribution-area with the greater Kye, 
and of Constituent Integral -bock @) from the cvertaping distribution-area — the lesser Kfe . 


pt a bi ‘ 
(2) When Prediction-area = 1,0000 (@) fep~eKy=4.71 R-units « T7151 - Gye 4128 (A560 - aa = .0356 = Constituent Integral- block 


fep~ cK,= 361 R-untts = 9246. O.,—»[.5872- (3224 -1.1744)] = .2086 = Constituent Integral- block 
(b) When Actual - arga = 1.0000 (3)fca~ Kg+9.71 R-untts « 1.4600: Og————» (4279 . 1.3333)=. 5705 
(.0828 sianie sood * AROL - Constrivent Integrat-biock @) 





fea~ cK, 1.39 R-untts = .2091-O%q 
p 2 


3 Constituent Integral- bioks ———> 1043= COEFFICIENT OF 
PREDICTION - ACCURACY (CPA) 
© LeGien equimelent number of Standard Deviations. 1. Numerical equivalent of area to left ~right,as the case moy be, o the given fc. 
I =Portion of total symmetrical area subtending given number of Standard Deviations. (From Table) 
WeTetal Area, if symmetrical,on basis of particular given J. 
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In this case the coefficient of prediction-accuracy is equal to the summa- 
tion of four constituent integral-blocks as follows: 


From that distribution-area with the greater K,,. 
Block 1 = area between « left and cKl. 
Block 2 = area between CKr and « right. 

From that distribution-area with the lesser K,,. 
Block 3 = area between cK/ and fc. 

Block 4 = area between fc and cKr. 


Case II is characterized by the set-up in which there are two common 
probability-values (cK/ and cKr). On the R-scale one such value lies either 
above or below and the second lies between the two fluctuation-centers 


(fep and fca). 


Kfe » 1.000 

—— 6065 
K On normal 
nN 


i<—K = .6065 


Normal Gaussian 


<—K<.6065 
Non normal 











Kfe = .0000 





>R fe 
FIGURE 3 
The effect, on ‘‘normality”’ and area, of making R = 2c, but in causing such R to cross 
the K-scale at any point other than at 0.6065 Kye. 





In this case the coefficient of prediction-accuracy is equal to the summa- 
tion of three integral-blocks as follows: 


From that distribution-area with the greater K,,. 
Block 1 = area between cK/ and « left. 
Block 2 = area between cKr and « right. 

From that distribution-area with the lesser K,,. 
Block 3 = area between cK/ and cKr. 


Case III is characterized by the set-up in which both the predicted and 
the actual areas possess the same-valued K,,, and each area is either bi- 
symmetrical or skewed, but if skewed the relationship between areas ‘‘p”’ 
and “‘a’’ is a right-and-left-mirrored skewness, that is, the right-hand stand- 
ard-deviation of one distribution-area is equal to the left-hand standard 
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CASE I. 
with Example. 


PROPERTIES (The Ay Prediction-area and the total Actual-area each =10000. Each is plotted 


the same coordinate frame ~ the same \K, that &, the same Freguency or Pro- 
babilitirsoales and the same R-class measure of the common subject 
1 R-class-range = 1 R-unit. { Area-unit = .01 K by 1R. 

toi Jype ty ht area-bounding curve, K= Kfc€ ex Fach curve with constant K fe eguad 
the other curve, and with two O% either egual to or different from each 

a but p Bea (G+) +2=O0 symmetrical. [If skewed, the two areas must be 
skewed in mirror- relation to each other, as in case of the given i 
(c) One common Probability-value, due to overlapping of curves "p" and "a". Jhis one cK 
ts always located midway between the two fc's ie. between feb and fea. 


(@) yf the particular probability curce, whether symmetrical or skewed, 
Are left fc: Area’ to right fc =: 0,. When total area = 1.0000, 
Area fo let fe = 4% Area to right fc» 


(e) CPA = Swo constituent Integral-blocks, one for oath probability-distribution-area. 
® i * for area with fc to right and @ £.,. 0“ for area with fc to left. 


GIVEN Prediction-Area = 1.0000 GIVEN Actual-Area - 1.0000 










K-SCALE 


FIND | 





at point 15.00 on R-scale ca at point 20.00 on R-scale 


s 0.0875 on K- scale - = 0.0875 on K-scalé 
~ 4.5593 R-units 2 4.5593 R-units 
= 4.2746 R- units ge = 48441 R-units 
s od 2 4.2746 R-units 


—————e 
<a 
ona 
K = Probability fe = fluctuation center CG = standard deviation symmetrical cKe >the left-most probability value 
R Measured Subject-quality 2 «= left of perticuler fe ~ 12. when Probability curve is common to both p and a 
p*e property of Prediction-eres 4 © right of particular fe ote. wel Seoed cK, *the corresponding right - most 
@ + @ property of Actual-area Kfc * probability value at particular fe Oe + OG = 205 ,when woth partic- Common’ probability 


ular area and Kfe ore constant. 9 * (ep ~fea) in R-units 


fe predicted ~ fc actual «g + 5.00 
R-iocus of cK when (fca~cK) © Xe, and (fep ~cK)=X, , Solve for the positive value of X, by means of the simultaneous eguations: 


Krep -€ > * Kfca €- The two formulae for the Common Probabtitty at cK, ai 

Xa9g-X~ Locus cK on R-scale « locus fe predicted +X,. X,*250 R-units Xgs - R-untts. 
The locus cK = 19.00 + 2.50 +17 50 
When the total area under the particular probability-curve = 1.0000 


(a) to left of fcp = 20ip+(M +O.) = .93754 +2 = 4688 1.0000 
(bo tb right of fcp = ode +04 oh) -1.00nte 0 - O58) : 

(c) ce of fca = 20%q (Oz +0.) Se. ee 1.0000 
@ toright of fca = 20,a—(&+Q) = .93754-+ 2 = 4686 


The of the equal Constituent Integral-blocks @) and ® from the overlapping areus of “p” and *a”, 

wll | ° I Iv 
(a) When Prediction-areg = 1.0000 fep ~cK +250 R-units = 5161 Op —>[5312- (1971 -10625)] = 3218 = Constituent Integratl-block 
(b) When Actual ~ area « 1.0000 fea~cK=250 R-units * 5161 -Ozq —> [5312 - (1971 - 1 0625)] » 3218 + Constituent Jntegral-block 


£2 Constituent Integral-biochs ————> 6436+ COEFFICIENT OF 
PREDICTION- ACCURACY (CPA) 
= Given equivalent number of Standard Deviations. 
= Numerical equivalent # Ara to the left or r right, as the case may be, of the given fc. 





L 
x 
I = Portion of total sy trical area ig given number of Standard Deviations (from Table). 
V =Total Area, if symmetrical, on basis of particular given O 
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deviation of the second, and vice versa. In this case there is only one 
common probability-value and this is always midway between the two 
fluctuation-centers (fcp and fca). 

In this case the coefficient of prediction-accuracy is equal to the summa- 
tion of two integral blocks—one from each of the two basic distribution- 
areas. 

Block 1, from the distribution-area with the K,, to the right, equals area 
between cK and « left. 

Block 2, from the distribution-area with the K,, to the left, equals area 
between cK and « right. 

The solution of these problems is made possible by use of the following 
properties of the probability-distribution area, that is, of the area under 
the Gaussian or normal curve of error: 

(1) The first basic probability-formula for the symmetrical probability- 
area, in which K = class probability and R = class-range in subject-value, 
is 

R?2 
K = f(R) = Kye ?**. 

The value of the probability K may be continuous, but R, the measured 
quality under consideration, must always be broken into a series of classes 
or grades. 

g = the standard deviation, e = 2.7182 which is the base of the system of 
natural logarithms. fc is the fluctuation-center of R-values. 

(2) The second basic formula for the symmetrical probability-distribu- 
tion area represents the value of the maximum ordinate as 

N 
Kp —- 

When distribution-data are given it is always possible to reduce the ab- 
solute number within each R-class to the percentage basis, and thus JN, or 
the total area, is made to equal 1.000. 

+/ 2 = 2.5066 is a constant, and when N is constant, K,, and o are in- 
‘terchangeable and are inversely proportional to each other. If Ky, be 
constant JN or total area varies directly with o, or if « be constant N varies 
directly with K,. (See Fig. 1.) 

(3) Regardless of the skewness of the particular distribution, graphi- 
cally the standard deviation crosses the K, at a K-point always equal to 


0.6065 times K,,. This value is deduced as follows: 
R2 
Given K = Kye 2% 
Let Ky, = 1.000 and let R = o. 


R2 
Then K = le 2? = z = a = 0.6065 which is a constant. (See 
ef Ve 
Fig. 2.) 
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(4) Further significance of the point 0.6065K on the maximum ordinate. 

This point 0.6065K,, is the only point on the maximum ordinate at 
which, if the K,, remain constant, the doubling of the R-value (in this 
case the a) will result in the exact doubling of the newly plotted distribu- 
tion-area. Thus, if R equal two original standard deviations, and if the 
point of maximum slope on the boundary curve be above 0.6065, the result- 
ing ‘‘non-normal”’ area will be more than doubled; if such point be below 
0.6065 the resulting area will be less than doubled—compared with the 
original normal Gaussian area for which the standard deviation equals 
1.000. (See Fig. 3.) 

Again, if the R-value of the original normal area be selected at less than 
one standard deviation (graphically above 0.6065K,.), and such R-value 
be doubled but kept at the same K-point (above 0.6065K;,,), the change 
causes a new ‘“‘non-normal”’ slope in the new boundary-curve-of-distribu- 
tion, and the new area will always be less than the original normal area in 
which the standard deviation equals 1.000. 

Or if the new R-value be greater than one standard deviation (graphi- 
cally below 0.6065K,.) and such R-value be doubled, and the new ‘‘non- 
normal’ area be plotted in the same manner as when the new standard de- 
viation is less than one, such new area will always be greater than the 
original normal area in which the standard deviation equals 1.000. (See 
Fig. 4.) 

Thus when the R-value = lo, and the K value = 0.6065, all other fac- 
tors remaining constant, such R-value is the only R-value, the doubling or 
other increase of which will cause a proportional increase in the normal dis- 
tribution-area. 

Next to the rectangle, the ‘‘normal-probability-distribution-area” is 
probably the simplest flat geometric figure the area of which varies directly 
with the change (that is, with the first power) of either one of its two co- 
ordinate factors, while the second remains constant. But there is this dif- 
ference: Plot any rectangle and any normal-probability-distribution-area 
each within an x-y frame of coérdinates, the area of the rectangle varies 
directly with the value of x selected at any y-point, or with the value of y 
selected at any x-point; while the area of the normal-probability-distribu- 
tion-figure varies directly with x only when x equals the standard deviation 
(which always crosses the y-coérdinate at the point equal to 0.6065 of the 
maximum ordinate), and directly with y only when y is the maximum ordi- 
nate, which always crosses the x-coérdinate at the fluctuation-center of x- 
values. 

(5) Relative areas. If the total normal distribution-area be 1.000, and 
the K,, be constant at 1K, then that portion of the area to the right 
of the particular K,, is to that portion of the area to the left of said K,, as 
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the particular standard deviation right is to the particular standard devia- 
tion left. 

Area right : area left = o, : oy. 

(6) Skewness. Using the second basic formula above given, we find 
that, if the K,, be constant the area varies directly with the standard devia- 
tion in a symmetrical or non-skewed area. Now, if the area be skewed, a 
hypothetical area must be postulated temporarily. Let, for example, the 
“standard deviation symmetrical” (o,) equal 0.50, then shift the “‘standard 
deviation left” (o,) to the value of 0.55. | Now if the K,, be constant at 
1K, and the area be made symmetrical with the new o = 0.55, the new 
total postulated-area equals X in the ratio 1:0.50 = X:0.55, in which case 
X = 1.10, ie., the oj would belong to a total area equal to 1.10 with its 
o symmetrical at 0.55 and its K,, constant at 1K—the original value. 


Kfc = 1.0000 


<—K > .6065 
Non- normal 





—K = 6065 
Normal Gaussian 


—K<.6065 
Non- normal 











Kfc = .0000 





fe 





FIGURE 4 


The effect, on “normality” and area, of doubling the normal R-value at any point 
on the K-scale other than at 0.6065 Kye. 


But actually the total area must remain constant at 1.000, and Ky, must 
remain constant at 1K. Then since the actual-skewed-area to the left 
has already consumed 0.55 of the total-actual-area, this leaves 0.45 for the 
area to the right of K,,. Similarly computed the o,;,n4 would belong to a 
newly postulated symmetrical distribution represented by a total area 
equal not to 1.000 but to 0.90, and with K,, constant at the same 1K 
value. 

Therefore 20 symmetrical = Otete + Tright- 

(7) Down the long S-curve (one of which bounds the top side of the 
probability-distribution to the left and another to the right of the particular 
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K,,), beginning at the top, where the slope equals zero, and descending the 
slope increases continuously until the descent-point, located at the height 
of K at lo removed from the particular fluctuation-center, is reached, i.e., 
until the descent-value 0.6065K,, is reached. This descent-point locates 
also the maximum slope of the long S-bounding-curve. (The greatness of 
this maximum slope is a function of its adjacent standard deviation c.) 
Therebelow the rapidity of descent decreases continuously, and for the 
latter reason never reaches the level whereat K = 0, but theoretically 
extends indefinitely, continuously approaching but never attaining 
parallelism with the particular R-scale at K = 0. 

(8) The total area under each curve (the prediction-area ‘‘p’”’ and the 
actual-area ‘‘a’”’) independently is equal to 1.000. The total area, including 
the overlap covered by both curves, equals (2 — CPA or overlap). 

The prediction-basis with the specifications listed in the earlier 
paragraphs is the one actually used. The selection of the particular pre- 
diction-basis must, of course, be justified by the logic of the particular case 
under analysis, but how the particular prediction-basis was derived or why 
used are matters which do not concern the basic theory of the coefficient of 
prediction-accuracy nor the technique of its use. The prediction-basis 
may thus be a guess or an intuition; it may be based upon a single proba- 
bility-distribution composed of a few individual cases carelessly measured, 
or it may be based upon the probability-resultant of many constituent 
probability-factors, each based upon many individual cases, and each ac- 
curately and confidently measured. An example of the latter type of pre- 
diction-basis is one produced by the probability-resultant* made by the 
proper synthesis of one cross-section from each of a number of selected 
manerkons each of which individually pictures mathematically the rela- 
tionship K = f(M, R), in which M is the particular prediction-basis, R is 
the thing-predicted and K is the probability. The point is that, in this 
particular study, the prediction-basis is being tested for its accuracy, not 
for its source. 

But in interpreting the meaning of a computed value for the coefficient 
of prediction-accuracy, the situation is such that one can criticize usefully 
the soundness of construction of the particular prediction-basis. It is evi- 
dent that if for many quantitative distribution-data on the same trait or 
quality, within the same population, the prediction-accuracy in many tests 
runs consistently high, there must be something besides luck or chance 
which favors the particular prediction-set-up. If the high running coeffi- 
cient of prediction-accuracy be based upon many independent factors, and 
each, in turn, is composed of many carefully measured individual cases, 
then confidence is justified that the essential major factors are in hand and 
are properly weighted and coérdinated in the given prediction-basis. 

But if the coefficient of prediction-accuracy run consistently low then, 








7 _ we 
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regardless of how carefully the particular prediction-basis is computed, 
it indicates that some major factor or factors are omitted, or are 
improperly weighted or coérdinated, in the perticular prediction-basis being 
tested for accuracy. 

In comparing two different prediction-bases against the same subsequently 
determined set of actual distribution-values, the coefficient of prediction- 
accuracy tests the relative merit of these two prediction-bases as to their 
relative correctness in the inclusion of their own postulated constituent 
factors, and in the relative weights and coérdinations given such factors. 

In any case in which one has in hand a prediction-basis in terms of a 
probability-distribution, and has also in hand the actual subsequent or 
other independently secured measuring and counting of the same quality 
in the same population, the coefficient of prediction-accuracy (CPA) 
presents a sound formula for judging the given probability-distribution- 
production-basis as to its degree of success—no excuses taken. 

The next step in good research is, of course, to find out why the particular 
prediction-basis which was actually used was not more successful. 


* H. H. Laughlin, “The Probability-resultant,”’ Proc. Nat. Acad. Sci., 21, 11, 601-610 
(1935). 


STRENGTH-DURATION CURVES AND THE THEORY OF 
ELECTRICAL EXCITATION 


By W. J. Crozier 
BIOLOGICAL LABORATORIES, HARVARD UNIVERSITY 


Communicated December 22, 1936 


I. An instructive parallel exists between measurements of electrical ex- 
citability, particularly for repetitive condenser discharges or constant cur- 
rent pulses, and the data of visual flicker. The comparison derives from 
probability considerations, and provides a means of analysis independent 
of any specific assumptions as to the physico-chemical character of the 
excitation process. 

Numerous observations describe the general relationship between mag- 
nitude of stimulating potential or current and the time for threshold excita- 
tion, particularly in nerve. It is well known that for a chosen degree of 
excitation the required current (or voltage) declines to an asymptotic level 
as the exposure time is increased.! No fully satisfactory theory of the 
form of the strength-duration curve has been provided. From a rational 
conception of the effects connected with electrical excitation Hill* derived 
equations describing rather closely the relation between current or poten- 
tial (C) and its duration (¢) for constant response. It is not clear that this 
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equation is an unequivocal fit for all series of observations. Its test in- 
volves visual adjustment of a curve on a log-log grid, since the equation 
cannot be conveniently rectified. There is some indication that the de- 
viations apparent in various sets of measurements, though slight, may be 
systematic. 
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FIGURE 1 
The relation of critical illumination for extinction of 
visual flicker (f) or for threshold response to flicker (F) is 
adequately described by a logarithmic probability integral. 
Only illustrative measurements involving cone excitation 
are here considered; for rod excitation a similar graph ob- 
tains.*? Ordinaic, 100 F/F,,,,; abscissa for A, log I 
(millilamberts), for B, log P (photons). Beyond J,,,,, F 
tends to fall, but the reasons for this® do not interfere with 
the present discussion. 
A. Data on the fish Xiphophorus; F,,,, = 43.1/sec.; 
21.5°C.; averages of 30 measurements. 
B. Three sets of measurements on the human eye (S. H.) 
from Hecht, Shlaer and Smith:* 
a—0.3° spot, white, on fovea (faz = 54.4) 
b—19° spot, 490 mu (fing, = 64.7). 
c—19° spot, 670myu (fmoz = 67.5). 
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Regardless of the descriptive adequacy of Hill’s equations, however, the 
data may be viewed from another standpoint. The general form of the 
relation between Jog C and log t is like that obtained in studies of visual 
excitation.* Here intensity (J) is the variable corresponding to ¢, and some 
measure of effect E corresponds to C. For threshold recognition of visual 
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flicker the critical flash frequency F is a proportionate measure of the in- 
tensity of the relevant total photosensory effect obtained.‘ Time is then 
involved only in the sense that the differential effect producing the re- 
sponse must be achieved (at any fixed flicker-frequency) within the interval 
set by the flicker cycle. In a voltage duration experiment a fixed level of 
effect is selected (threshold response, or half-maximum of action current, 
for example), and the stimulating current C (of fixed frequency) requires a 
certain measured time ¢ to produce this effect. The two kinds of measure- 
ments are therefore alike in that of the three factors involved—time of 
action, intensity and level of total effect—one is in each case essentially 
fixed by the procedure. 

The stimulation effect achieved is in each case due to the cumulative 
action of a number of excited elements—in a nerve or muscle either a num- 
ber of different fibres or a number of “‘excitable spots’’ on each fibre or both; 
and in the reaction visually excited, a number of retinally connected units 
must be excited. These dynamical units or elements are not to be identi- 
fied with the crudely obvious structural units, nor are they to be conceived 
as morphologically fixed. The excitability of each element is not station- 
ary but fluctuates as a function of intensity and of time. When E is arbi- 
trarily chosen, the time required for its production will depend upon the 
time-summation of contributions from a number of elements forming part 
of a frequency distribution of time X intensity thresholds. In such situa- 
tions one expects the law for E as a function of J (or #) to be dependent 
upon a logarithmic probability distribution of excitability thresholds, the 
latter given by (¢)(d@E/dt). This is quite without reference to the intimate 
physical basis of the exciting process in any given element.® 

The curve for critical illumination as a function of flicker-frequency® 
is admirably described as a log-probability integral (Fig. 1). Equations of 
the photostationary state type,’ with certain limitations, do so equally 
well because they are identical with the logarithmic logistic, which except 
at the extreme ends of the curve gives a fit scarcely distinguishable® from 
that of the log-probability integral; however, the proportionality constant 
in the “‘photostationary state’’ equation has properties which forbid assign- 
ing this particular significance to the equation.® The log-logistic equation 
also describes the measurements of cochlear effect as a function of intensity 
of sound (pressure). 

II. Strength-duration curves for repetitive electrical excitation are very 
efficiently described by a log-probability integral. Several examples are 
given in figure 2. Series of measurements have been selected from the 
most recent work which seem to be intrinsically the best determined and 
which cover the largest range of values of C and ¢. Numerous less exten- 
sive sets of data!!? are equally well described by this equation. The test 
is easily made on a log-probability grid (Fig. 2). This requires that the 








74 PHYSIOLOGY: W. J. CROZIER Proc. N. A. S. 


acting current (or voltage) be computed as a percentage of the asymptotic 
minimum Cy. That is, 100 ()/C represents the percentage of available 





99.9 
99 


35 
90 
80 


79 


10 OfC 
° 
a “SS 














0.05 


fog ¢ 


FIGURE 2 


Log probability grid: ordinate, 100 C,)/C; abscissa, log RC, the condenser 
time constant. Repetitive condenser discharges. 

(i) Data from: 4 experiments in Hill, 1936);? sciatic of frog; II, 5 shocks/ 
sec., 4°C. (Co = 0.248). Log#’ = 3.35. IIIa, 31.4 shocks/sec., 22.5°C. 
(Co = 0.708). Log #’ = 4.84. IIIb, 17.2 shocks/sec., 6.2°C. (Co = 
0.347). Log t’ = 3.245. IV, 31 shocks/sec., 20°C. (Cy = 0.233). 
Log t’ = 3.04. The slope is the same for each line. 

(ii) Data in: Scott, D., 1934, Jour. Physiol., \xxxii, 321; frog sciatic, con- 
denser discharges 9.6/sec. (Co = 0.52). Log #’ = 2.79. [Data read 
from curve.] 

(iii) Data from: Rosenberg, H., 1935, Jour. Physiol., \xxxiv, 50. a—cat, 
sciatic, 21.2° (Cy = 0.2296). Log t’ = 2.74. b—crab (Maia), leg 
nerve, 3 mm. electrode separation, (Co = 0.378) and 15.5 mm. (Cy = 
0.320). Logit’ = 5.25, 5.40. 

(iv) Data from: Bugnard, L., and Hill, A. V., 1935, Jour. Physiol., Ixxxiii, 

' 4265, fin nerve of Sepia. 
a—I, 1’ electrode separation 1.7 mm. (Cy) = 5.76). Log t’ = 3.285. 
I, 1” electrode separation 16 + 3 mm. (Cy = 5.22). Log t’ = 3.44. 
b—I, 2’ electrode separation 1.7 mm. (Cy = 3.17). Log t’ = 2.86. 
I, 2” electrode separation 13-19 mm. (Cy = 2.7). Log t’ = 2.93. 

(v) Data from: Jinnaka, S., and Azuma, R., 1922, Proc. Roy. Soc. (London) 
B, xc, 49; frog sartorius muscle; scale of log ¢ at the top. 
a—Single fibre, pore electrode (Cy = 1.063). Log é’ = 4.519. 
b—Pelvic end of sartorius, with 0.05 per cent curare (Cp = 0.730). Log 

t’ = 4,82. 


“elements,” in terms of the intensity X time relationship, stimulated by C 
at the moment when threshold response occurs. Put in another way, for 
very short times 1/C measures directly the excitability of any element of ex- 
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citation, and should give a normal frequency distribution of excitabilities. 
At the lower limit of 1/C this gives the “‘constant quantity” rule. Cy is, 
of course, not to be confused with the observed rheobase, but is a little 
below it and must usually be found by trial. The graph is quite 
sensitive. 

This treatment discloses that the voltage duration data in experiments of 
one type with one kind of nerve often may be compared in terms of a single 
parameter, namely, the value of log ¢ for which Co/C = 50 per cent (cf. 
Fig. 2). For example, in four series of measurements from Hill’s recent 
paper! the values of RC = ¢ have been multiplied by a constant for each 
series; the four sets are seen (Fig. 3) to be identical as regards the shape of 
the curve. The value of ¢ for Co/C = 50 per cent is a function of tempera- 
ture, of interelectrode distance, of fibre diameter and of the individual 
preparation. With fall of temperature Co/C = 0.50 moves to a higher 
value of ¢, as is expected on the basis of a frequency distribution of excita- 
bilities. The effect of electrode separation,? although probably more com- 
plex, is also of this general character; C falls, t for C = 2Cy rises and o,, , 
is increased, as it should be if a larger population of excitable elements is 
made available. The parallels with the properties of flicker data are pre- 
cise. For different types of nerve and muscle and for other differing ex- 
perimental conditions, the standard deviation of the underlying frequency 
distribution of excitabilities is not the same, so that an additional parame- 
ter (Gig ,) is needed. Its properties should repay systematic investiga- 
tion. Graphs with different slopes might be compared by finding that ad- 
justment of the log ¢ scale units which will give them the same slope, after 
rotation about the 50 per cent point, by employing log # instead of log ¢ 
and having ” proportional to the slope on the probability grid. 

III. The value of ¢ for 100 C)/C = 50, the inflection point of the curve, 
has a special significance. It gives the time for excitation by a current = 
2Co, under given conditions, and consequently corresponds in meaning to 
the chronaxie of Lapicque'*—although computed on the basis of the 
asymptotic level of Cy and not from the observed apparent rheobase. 
The curves show that this ‘‘true” chronaxie, which may be called 1’, is a 
valid parameter of the excitability function. It is independent of Cp. 
It is not dependent upon the position of the curve at one point, but (cf. 
Fig. 2) is weighted by all the observations. The physical meaning of r’ 
is the time required for excitation by that current which under the condi- 
tions imposed stimulates up to the median (mode) of the frequency distri- 
bution of excitability thresholds. This relationship between r’ and the 
classical chronaxie or excitation time is curious because it was largely a 
matter of accident that twice the rheobase voltage should have been chosen 
originally for the measurement of the excitability time factor.'® 

It is obvious that 7’ alone is not adequate for the efficient description of 
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FIGURE 3 
The four curves of (i) in Fig. 2 are shown brought together at C)/C = 50 per cent, RC 
being multiplied by the factors 1, 3.273, 1.35 and 2.19, respectively, for II, IIIa, IIIb 
and IV. The curve is that for the probability integral 
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the strength-duration curve.® The slope factor (i.e., o,, ,) is also required 
if different kinds of preparations are to be compared. 

With these two parameters it is not difficult to compare different forms 
of excitability, for example, photic excitability (Fig. 1) and electrical (Figs. 
2, 4) in diverse irritable structures. In a similar way the data of visual 
excitation may be investigated as a function of wave-length with the expec- 
tation of evidence as to the existence of different kinds of receptor excita- 
bility distributions. The formulation proposed is in terms of the fre- 
quency distribution of excitabilities among a population of excitable ele- 
ments and is independent of assumptions as to the nature and the inter- 
relationships of the physico-chemical changes which produce excitation in 
any one kind of sensory element. The result is thus consistent with that 
arising in the analysis of the sensory discrimination of intensities. '” 

IV. Summary.—Measurements of the interrelations of current strength 
and time for a fixed magnitude of effect due to repetitive electrical excita- 
tion in nerve and muscle are found to be described with excellent fidelity 
by a logarithmic probability integral. The same kind of relationship, in 
which magnitude of sensory effect is expressed as a function of intensity, 
appears in phenomena of stimulation by light and by sound. The value of 
time for which 100 Co/C = 50 per cent, where Cy is the asymptotic lower ex- 
citing current, is essentially the chronaxie. This parameter is not by it- 
self sufficient for the description of the strength-duration curve; the stand- 
ard deviation of the underlying frequency distribution is an independent 
parameter, being different in various kinds of nerves and according to con- 
ditions of excitation. The log-probability integral gives a readily applied 
means of comparing different forms and manifestations of excitability. It 
is pointed out that, when an exciting energy is applied to a tissue, it is to 
be expected that the thresholds for excitation of the elements open to 
stimulation, in terms of the resultant chosen, will form a logarithmic fre- 
quency distribution provided the excitability of each element intrinsically 
fluctuates. 
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as they go, are quite well described by a log probability integral. 
In the iron wire model used by Lillie (Jour. Gen. Physiol.. 19. 109 (1935-36)) the “ex- 








78 MATHEMATICS: H. RADEMACHER Proc. N. A. S. 


citation’”’ of a number of “irritable spots’ on the wire, resulting in the activation of a 
propagated wave, will scarcely be influenced by fluctuations in catalytic activity at each 
spot. With constant current, the voltage duration data for activation are well de- 
scribed by a simple probability integral. 

6 Crozier, W. J., Wolf, E., and Zerrahn-Wolf, G., Jour. Gen. Physiol., 20 (in press) 
(1936-37 8, d). 

7 Hecht, S., Shlaer, S., and Smith, E. L., Cold Spring Harbor Symposia Quant. Biol. 
3, 237 (1935). 

* Cf. Wilson, E. B., Proc. Nat. Acad. Sci., 11 (1925); Winsor, C. P., Jour. Wash. 
Acad. Sci., 22, 73 (1932). 

® Crozier, W. J., Wolf, E., and Zerrahn-Wolf, G., Jour. Gen. Physiol., 20 (in press) 
(1936-37 a, b). 

10 Cf. Crozier, W. J., and Holway, A. H. (in preparation) (1937). 

11 Cf, 1} 2and Davis, H., and Forbes, A., Physiol. Rev., 16, 407 (1936) for citations. 

12 In many instances (e.g., nerve-muscle preparations and others) two intersecting 
strength-duration curves are obtained due to the presence of two groups of excitabilities. 
In some cases the separation may not be clear cut. The analogy with the rod and cone 
situation in the excitation of the vertebrate retina is obvious. 

13 Hill, A. V., Proc. Roy. Soc. (London) B, 119, 440 (1936 b). 

14 Lapicque, L., L’Excitabilité en Fonction du Temps, Paris, Presses Universit., 
371 pp. (1926). 

15 Lucas, K., Jour. Physiol., 35, 326 (1906); Lapicque, L., C. R. Soc. Biol., Paris, 67, 
280 (1909). 

16 Cf. Davis, H., and Forbes, A., Physiol. Rev., 16, 407 (1936). 

17 Crozier, W. J., Jour. Gen. Physiol, 19, 503 (1935-36); Proc. Nat. Acad. Sct., 22,412 
(1936); Upton, M., and Crozier, W. J., Ibid., 417 (1936); Crozier, W. J., and Holway, 
A. H., Ibid. (in press) (1937). 


A CONVERGENT SERIES FOR THE PARTITION FUNCTION 
p(n) 
By Hans RADEMACHER 
DEPARTMENT OF MATHEMATICS, UNIVERSITY OF PENNSYLVANIA 


Communicated January 9, 1937 


I. For a discussion of the number p(m) of unrestricted partitions of ” 
there has been used since Euler! the generating function 


oO 1 
fx) =1+ >» b(n)x* = (1 — x)(1 — x*)(1 — x)... (1) 





A purely formal treatment of (1) leads to classical formulae of recurrence. 
In 1917 Hardy and Ramanujan? applied on (1) the theory of functions of 
complex variables and developed a method (Cauchy integral, “Farey- 
dissection,” approximation functions) which yields an asymptotic formula 
for p(n). Their result is the following: 
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CVn — 'u 
1 ae . = 
p(n) a QavV/2 Zim Vie “mae + O(n 4) (2) 


with a as an arbitrary constant, and 


— 2rthn 
A,(n) = » * Wh, k & “ »,C=4rvV s/s. (3) 
mod k 


anes 


The symbol w, , designates certain 24kth roots of unity which have their 
origin in the theory of the transformation of modular functions.? The 
formula (2) is noteworthy for the fact that its error term tends to zero with 
increasing n. Tests for m = 100 and m = 200 showed indeed that p(n) 
can be obtained as the nearest integer to the sum expression in (2) for a 
small number of summands. 

The following revision and modification of the Hardy-Ramanujan 
method permits us to replace the asymptotic formula (2) by the equality 


q sinh S¥* _* 
p(n) = mre i A,(n)Vk > woe Pe (4) 


in which the series is absolutely convergent. The asymptotic formula (2) 
is then an immediate consequence of (4). 
By Cauchy’s theorem we deduce from (1) 


-iJ: LO) ae, (5) 


where C is a certain circle enclosing the origin and interior to the unit 
circle. The approach to a treatment of (5) is given by the formula 


Qxih’ Qn 


NaF) 0 wae (<5. - ays sl 


with 2(z) > 0; ~/z denoting the principal branch of the square root and 
h’ determined as a solution of the congruence hh’ = —1(modk). The 
formula (6) itself is a consequence of a transformation formula for 
n(r), a modular function introduced by Dedekind,‘ which is essentially the 
24th root of the “discriminant.” 

In (5) we choose C as the circle of radius e~**"* and divide it by the 
Farey-dissection of order N. The Farey arc & , may be described by 


anit — 2xN~-? + 2rig 
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Introducing (6) into the integrand we get 


1 2 f(x) 
njp=>—.,. h,k)= es 
b(n) Omiodte ny J xetids 
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where we have used the abbreviation 


eh a Sa (=. =) - 
P \idke 12k ” 
Instead of introducing a special auxiliary function in order to approximate 
the integral, as Hardy and Ramanujan do, we use V(z) itself, since f(x’) 
approaches 1 very rapidly when x’ tends to zero. If we carry out the 
proper estimates we obtain 





2Qxth dh, 
wenn 
iar oe ee 7 U(R(N ~* — ig)) e-™*™* dp + 0, 
hk =1 
OSh<kSN —9;, k 


acer" a. (8) 


It may be noticed that in our reasonings the order N of the Farey-dissec- 
tion is independent of the represented number . 

In (8) we substitute w = N-* — ig and consider the integral in the 
complex w-plane. In order to get a unique determination of V(kw) we 
cut the w-plane along the negative real axis. The integral in (8) is then 
replaced by a loop-integral extended over a path surrounding the cut and 
connecting the following points by straight segments: 


*of *of i +o! *o/ 
Sy ee Wh, ks aN 19h, ks itis + Wh, ky ~~ > Dh, ks as ee 


The estimate of the errors is based on the fact that the essential singularity 


of ¥(kw) in w = 0 arising from a factor e” with positive a has bounded 
values for R(w) < 0 and that more precisely for w = u + iv we have 


for u < b(u? + v?), which, for b > 0, is fulfilled a fortior: if u S bv?. This, 
by the way, is the reason for choosing the radius e~ **N™ together with the 
order N of the Farey-dissection. 
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The calculation leads to 


N 
p(n) = Y Aaln)ox(n) + Or (9) 


| Or| < Bren, 


where A,(m) is the quantity defined in (3) and 
__ (0+) 
ox(n) = “> - Vw elena 9" dry 


= @ 


+ VE f Ve BO ay, (10) 
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The remainder term in (9) tends to zero for N —>o and fixed nm. Hence 
we have 


p(n) = DALnor(n), (11) 


The integrals in (10) are well known, and we get in this way the proof of 
(4). 

From (4) it is possible to derive a definite expression for the O-term in 
(2). Weobtain 


CVn—Vu (12) 
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a formula which proves to be useful for numerical computation of p(m). 
A detailed account of the outlined proofs will soon be published in the 
‘Proceedings of the London Mathematical Society.” 


II. The result (4) can in turn be used for a treatment of f(x), which 
because of its close relation to modular functions deserves some attention 
for itself. We get from (1), (4) and (3) 
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is an entire transcendental function in of order 1/2. Hence by a theorem 
of Wigert® the function 


sinh 





u(e) = aa(ns" 


originally defined only in the interior of the unit circle, can be continued 
over the whole z-plane, is regular at z = ~ and has only one isolated essen- 
tial singularity at z = 1. 

The expansion (13) can therefore be taken as a generalized Mittag-Leffler 
decomposition into ‘‘partial fractions," each term giving rise to an essential 
singularity at a root of unity x = e*. The particularity of (13) lies 
in the fact that all singularities together build up the natural boundary 
|x | = 1 for the function f(x). 

The ®,(z) have, by Wigert’s theorem, also a meaning outside the unit 
circle, and a simple consideration shows that the right-hand side of (13) 
is again convergent for | x | > 1. This defines there an analytic function: 

2nih 


fe) =14 25S  onevin(xeF) lel > 1. (30) 


(h, k)=1 
h mod k 


The functions f(x) and f(x) are determined by the same expansion (13), 
(13a), in spite of the fact that they are separated by the unit circle as natu- 


ral boundary of f(x). 
We can establish a power series for f"(x) aroundx =. From (14) and 


(15) we get 
oe (¢ A 
o\ek p> © 
$,(z) = a (+1)! DIY (n — t/,)’~" s n => “ o> FiO) 
say. Wehave 
ax Sr = =— -- ds" 








ward —2Z 1—27} 








‘it 
3) 


3), 
tu- 


ind 





VoL. 23, 1937 MATHEMATICS: H. RADEMACHER 83 


and 
ots) = (64 — Yu) os) 
“ dz 
Consequently around z = ~ the expansions 
en(2) = (17 Do (m+ Yu) 2, || > 1 
are valid, and we have 
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with the trigonometric sine instead of the hyperbolic sine. With this 
determination (13a) takes the form 
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A formal comparison of (4) and (17) can be condensed in the equation 
p*(m) = —p(—m). 
Whereas p() is very rapidly increasing with increasing m, p*(m) turns out 
to be rather small. An estimate readily obtained is 
p*(m) = O(m~ ). 


The function f*(x) requires a further investigation. It is quite possible 
that f*(x) is identically zero for | x | > 1. Not only certain function-theo- 
retical considerations, but also some numerical computations pertaining 
to the first coefficients p*(m), which Dr. D. H. Lehmer was kind enough 
to carry out, seem to point in this direction. 
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The whole theory admits an application to a class of functions related to 
certain ‘“‘modular forms’’ of positive dimensions, e.g., to f’(x) for 2 S r S 24. 
Dr. Zuckerman and I are studying these generalizations in a paper now 
under preparation. 


1 Introductio in analysin infinitorum (1748), cap. 16, ‘“De partitione numerorum,” 
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tische Funktionen, 97-102 (1891), and H. Rademacher, Jour. reine angewandte Math., 
167, 312-336 (1931). 

5 “Sur un théoréme concernant les fonctions entiéres,’’ Arkiv Mat., Astron. Fystk, 
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NOTE ON THE CURVATURE OF LEVEL CURVES OF GREEN’S 
FUNCTION 


By J. L. WALSH 
DEPARTMENT OF MATHEMATICS, HARVARD UNIVERSITY 


Communicated December 19, 1936 


The primary object of this note is the proof of 

THEOREM 1. Let the function w = f(z) with {(0) = O map the closed 
region \2| < 1 smoothly onto the closed interior of a convex analytic Jordan 
curve C in the w-plane. Let C,,O <r S 1, denote the image in the w-plane 
of the circle \2| =r; that is to say, C, ts the level curve G(x, y) = log r of 
Green’s function G(x, y) for the interior of C with polein w = O. Then the 
expression 

r- [min. curvature of C,] (1) 


increases monotonically as r decreases, and approaches the limit 1/ lf ‘(0)| as fr 
apbroaches zero. 

For the present, the term ‘‘increases monotonically’’ is to be interpreted 
in the weak sense ‘‘increases or remains unchanged.’ 

Let 2 with l2| = r(O <r & 1) bearbitrary, and let x(z) denote the curva- 
ture of C, at the point w = f(z). We state for reference the well known and 


easily proved formula 
swf +st | 
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where the symbol R denotes the real part of the complex quantity which 
follows it. Equation (2) is to be interpreted algebraically rather than 
numerically; if z traces the circle |2| = rin the positive sense, then x(z) is 
positive or negative according as the tangent to C, at w = f(z) rotates in the 
positive or negative sense. 

Under the hypothesis of Theorem 1 the curvature of C is non-negative 
and continuous. If the curvature vanishes on C, it cannot vanish at every 
point of the Jordan curve C, and is positive where it is not zero. Then the 


harmonic function 
uw 
af 1 + of | 
Ff’) 
is actually positive at every interior point of r’;|2| = 1. The function 


f'(z) is continuous and different from zero on and within [. That is to 
say, if x(z) vanishes at points of I’, we have shown for \2| =r<l 





r-k(z) 2 min [k(z), zon I]. (3) 


It remains to prove (3) when «(z) does not vanish on I. 

We shall have occasion to apply the following lemma, whose proof is 
merely a matter of computation: 

Lemma. If the function u(x, y) 1s harmonic and positive, then for the 
function U(x, y) = log u(x, y) we have 


2 2 2 2 
OF ee ley +(24)']s0 
Ox? = Oy’? u?L\Ox oy 
Consequently! the function U(x, y) is superharmonic. 
We are now in a position to complete the proof of Theorem 1. Under the 
present hypothesis the curvature of C is positive, so the numerator in (2) 


is positive on I’, and being harmonic is positive also for lz| <1. Let us 
write (2) in the form 





log {lz|-x(@)] = log at| 142 | —loglf'@)|- (4) 
f'@) 

The first term on the right is superharmonic on and within I’, by the 
Lemma. The second term on the right is the negative of the logarithm 
of the modulus of a non-vanishing analytic function, hence harmonic on 
and within’. Then the entire right-hand member of (4) is the sum of two 
functions superharmonic on and within I, hence itself is superharmonic on 
and within. The function log [lz|-«(z)] is continuous on and within I, 
and being superharmonic can have no proper minimum interior tol. Then 
the minimum of log [z|-«(z)] or of |2|-x(z) in the closed interior of I is as- 
sumed on I itself, so we have established (3) for [z| =r<l. 
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The function w’ = f(rz’) with O <r <1 maps the region |2’| < 1 smoothly 
onto the closed interior of C, in the w’-plane, so from (3) we may write for 
le’] = n/r < lor |s| = 1 <r 


“x(e) 2 min [x(2), |2| = 1]; 


the monotonic character of (1) is now established. The remaining part of 
Theorem 1 is obvious by inspection of (2). 

Let K, be the curve obtained from C, (r S$ 1) by stretching the w-plane 
away from the origin O in the ratio r : 1. Theorem 1 asserts that the 
minimum curvature of K, increases monotonically as r decreases and ap- 
proaches the limit 1/ lf : (0)| as t approaches zero. Of course the maximum 
curvature of K, also approaches 1/|f(O)| as r approaches zero. It is also 
true whether C is convex or not, and is easily proved by studying the 
modulus of the function f(z)/z, that the greatest distance from O to a point 
of K, decreases monotonically as r decreases and approaches 1/ | f’ (0)| asr 
approaches zero; by a similar study of the modulus of 2/f(z) it may be 
proved that the least distance from O to a point of K, increases monotoni- 
cally as r decreases, and approaches 1/|f’(O)| as r approaches zero. 

We remark that the maximum curvature of K, does not necessarily de- 
crease monotonically as r decreases. Indeed, let C be the circle of radius 
unity whose center is the point w = a, O < a < 1, with 
s—a 


w= f(z) = 





as=— 


we easily compute for |z| = r the value 


— gr2 
nscale wo 


which depends only on r = {2|. Then the curvature of K,, which is con- 
stant on K, and equal to the maximum curvature of K,, increases mono- 
tonically in the strong sense as r decreases. 

We shall show later that for arbitrary convex C the maximum curvature 
of K, does not necessarily increase monotonically as r decreases. 

If the curvature of C is not actually positive at every point, the strong 
inequality is valid in (3). Suppose now that the curvature of C is positive 
at every point; we shall prove that except in a trivial case the strong in- 
equality is also valid in (3). 

The function log {lz|-«(z)] is continuous and superharmonic on and 
within T. Let V(x, y) denote the function which coincides with log 
[z|-«(z)] on I’, and is harmonic interior to [ and continuous in the corre- 
sponding closed region; the minimum value of V(x, y) on and within T 
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can be taken on by V(x, y) interior to T only if V(x, y) is identically con- 
stant. Thus the minimum value of the superharmonic function log 
{|z|-«(z)] can be taken on by log [z|-«(z)] interior to T only if V(x, y) is 
identically constant. In the latter case x(z) is constant on I, so C is a 
circle; the reader will verify from the actual formulas involved (see above) 
that the point w = O must lie at the center of C. Then in every case 
where C is not a circle whose center is w = O, the strong inequality is 
valid in (3), so under the hypothesis of Theorem 1, except in the trivial 
case that f(z)/z is identically constant, the expression (1) increases in the 
strong monotonic sense as r decreases. 

Under the conditions of Theorem 1, the minimum curvature of C is posi- 
tive or zero. If Cisa circle, each curve C,(r < 1) isconvex. If Cisnota 
circle, the strongly monotonic variation of (1) implies that the curvature of 
C, (r < 1) is positive, so C, is also convex. he conclusion that in every 
case C, (r < 1) is convex is a well known theorem due to Study,? thus con- 
tained in Theorem 1. 

Another generalization of Study’s result is the theorem proved by 
Peschl,* that if the function w = f(z) with f(0) = O maps |z| < 1 onto the 
closed interior of an analytic Jordan curve C whose minimum curvature is 
k, then the minimum curvature of C, is greater than or equal tok. Peschl’s 
result is more general than Theorem 1 in admitting either positive or nega- 
tive curvature of C, but is less general in that Theorem 1 asserts for k = O 
the monotonic character of (1). From the geometric viewpoint, Peschl’s 
result thus does not involve the monotonic character of the minimum curva- 
ture of K,. 

It is a consequence of Theorem 1, together with the additional remark on 
the strongly monotonic variation of (1), that under the hypothesis of Theo- 
rem 1 the minimum curvature of C, increases in the strong sense as r de- 
creases. It cannot be proved that the maximum curvature of C, always 
decreases as 7 decreases; indeed in the case f(z) = z the maximum curva- 
ture of C, increases as r decreases. 

It also cannot be proved that under the conditions of Theorem 1 the 
maximum curvature of C, always increases as r decreases. Let w = F(z) 
map |z| < 1 onto the interior of a lens-shaped region D in the w-plane, sym- 
metric in the axis of reals and in the origin, bounded by two equal arcs of 
circles which intersect in points A and B. Let D, denote the image in the 
w-plane of the circle |z| = R< 1. The Jordan curves D, are convex curves 
whose curvature becomes infinite as w approaches A. Let R; be arbitrary, 
O < R, < 1, and let the maximum curvature of D,, be k;. Let R; be 
chosen, Re < R; < 1, so that the maximum curvature of Dg, is greater than 
ky. Then the function w = F(R’) maps |s’| < 1 onto a region bounded by 
the convex analytic Jordan curve D,,. Under this map the maximum 
curvature of D,,, the image of |2'| = 1, is greater than the maximum curva- 
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ture of D,,, the image of |2’| = R,/R, < 1. Wecan also choose the value 
of Ri, Re < Ri < 1, so that we have 


R,:[max. curvature of D,z,] > Re:[max. curvature of D,,], 


which shows that the maximum curvature of K, does not necessarily in- 
crease monotonically as r decreases. 

Theorem 1 is formulated for the case that C is analytic, but obviously 
applies in more general cases, for instatice whenever (2) is valid and x(z) is 
continuous for |z| < 1. Whenever «(z) is positive interior to I’, that is to 
say, whenever C is convex, we can replace (3) by the inequality 


r-x(z) 2 lim «(z), 
s—>T 





where the right-hand member represents the inferior limit of «(z) for all 
possible sequences of points z interior to T approaching I. 

If in Theorem 1 we omit the requirement that C shall be convex, and if 
curvature is interpreted numerically rather than algebraically, it is still 
true that (1) increases monotonically in the weak sense as 7 decreases. If 
C is not convex, expression (1) is then zero for r = 1 and remains zero for 
values of r for which C, is not convex. If 7 is the greatest r for which C, 
is convex (such a value 7; > O exists), expression (1) increases strongly 
monotonically as r S 7 decreases. 

Study of the function log [—x«(z)/le|], where x(z) is defined by (2), yields 

THEOREM 2. Let the function 


w= fie) = + at+az+az?+....0<0, 
Z 


map the closed region \z| < 1 smoothly onto the closed exterior of a convex 
analytic Jordan curve C in the w-plane. Let C,,O <r S 1, denote the image 
in the w-plane of the circle |z| =r; thatts to say, C, 1s the level curve G(x, y) = 
log r of Green’s function G(x, y) for the exterior of C with pole at infinity. 
Then the expression 





min. curvature of C, 
: (5) 


increases monotonically as r decreases, when curvature 1s interpreted numert- 
cally rather than algebraically ;> the increase 1s strongly monotonic if C is not a 
circle. Both (5) and 
max. curvature of C, 
r 





approach the limit 1/ la| as r approaches zero. 
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Let K, denote the curve obiained from C, by shrinking the latter toward 
the point® w = ao in the ratio 1:r. Then the minimum curvature of K, in- 
creases monotonically as r decreases; both maximum and minimum curva- 
ture of K, approach the limit 1/ la| as rt approaches zero. Whether C is con- 
vex or not, the greatest distance from ao to a point of K, decreases monotonically 
as r decreases, and approaches 1/ la| as t approaches zero; the least distance 
from ay to a point of K, increases monotonically as r decreases, provided ag 
is interior to K,, and approaches 1/ la| as t approaches zero. 

1 For this and the other elementary properties of superharmonic functions that we use, 
the reader may refer to Kellogg, Potential Theory (Berlin, 1929), p. 361. 

2 Konforme Abbildung einfach-zusammenhdngender Bereiche (Leipzig, 1913). 

3 Peschl, E., Math. Annalen, 106, 574-594 (1932). 

‘ But when 7 is sufficiently small, whether C is convex or not, the maximum curvature 
of C, increases as r decreases. 

5 But the minimum curvature of C, need not increase monotonically nor decrease 
monotonically as 7 decreases. 


6 The point w = do is the conformal center of gravity of the curve C, hence lies interior 
to C. See Frank, P., and Léwner, K., Math. Zett., 3, 78-86 (1919). 


REPRESENTATIONS INDUCED IN AN INVARIANT SUBGROUP 
By A. H. CLiFrForD 
INSTITUTE FOR ADVANCED Stupy 


Communicated December 21, 1936 


If h — A(h) is any representation %,, of a group H in any field K, A(h) 
indicating the matrix corresponding to the element h of H, and if h > h 
is any automorphism of H, then the representation 8, of H defined by 
h — B(h) = A(h) will be called an automorph of U,. If H is contained 
in a group G, and h = u—hu with u in G (H invariant under x), we shall 
say that %, is an automorph of ,, relative to G. 

If g — A(g) is any representation %, of a group G, and if H is any 
subgroup of G, then %, induces the representation Y%,, of H defined by 
h — A(h). Regarding the representation Y,, induced by an irreducible 
representation %, of G in an invariant subgroup H of G, we have the 
following general result: 

THEOREM 1. Y,, 1s either itself irreducible or is fully reducible into irre- 
ducible components all of the same degree, and all are automorphs relative to 
G of any one of them. 

Let Uy, Ay, . . .. Ay be a complete list of the inequivalent irre- 
ducible components of %,. Choose a codrdinate system in the representa- 
tion space ® of %, which exhibits the full reduction of %,, and let ®, 
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denote the sum of all the equivalent irreducible invariant subspaces of 
belonging to %,?, (i = 1, 2,..., 7). 

THEOREM 2. The spaces Ri, Re, . . ., R, constitute a system of imprimi- 
tivity of Ug. 

Since %, is irreducible, it is transitive, and all the spaces R, have the 
same dimension. Since all the irreducible components of &, have the 
same degree #, each of the inequivalent irreducible components of A, 
must occur the same number of times s, and the degree of %, is therefore 
rst. In particular, if U%, contains at least two inequivalent irreducible 
components (r > 1), then %, is imprimitive.! 

In the particular case when the factor-group G/H is a finite cyclic 
group, and the number field K is algebraically closed, we can give a rela- 
tively complete description of the situation. 

Let m be the index of H in G, and let us write 


G=H+wuH+wH+...+ vu", (uv in H). 


Let € be any mth root of unity in K. If Wg is any representation g — A(g) 
of G, and we define 

A*(uth) = &A(uth), (4 = 0,1,...,m — 1; hin A), 
then g — A*(g) is a representation %,* of G. We call %,* an associate 
of &,, and more specifically an associate of index k if ¢ is a primitive kth 
root of unity. 

With the aid of this definition, we may describe the situation as follows. 
All the irreducible components of %{, are inequivalent (s = 1), and their 
number 7 is a divisor of the index m of Hin G. %g is equivalent to all of 
its associates of index r, but to none of higher index (this serves to deter- 
mine 7; in particular, %,, remains irreducible if and only if %, is equiva- 
lent to none of its associates other than itself). Tf 6, is any other irre- 
ducible representation of G, and if one of the irreducible components of 
%, is equivalent to one of those of %,, then By is equivalent to A%,, and 
B, is equivalent to an associate of %,; the converse is evident. Finally, 
any irreducible representation of H can be obtained from an irreducible 
representation %, of G as one of the irreducible components of %,. 

These results were known by Brauer for the case when G is the full 
orthogonal group and H the proper orthogonal group,? and by Professor 
Hermann Weyl—to whom I owe the inspiration for the present develop- 
ment—for the case when H is a subgroup of index two in an arbitrary 
group G. Proofs will be given in a paper to appear in the Annals of 
Mathematics. 


1A. Speiser, Theorie der Gruppen von endlicher Ordnung (Berlin, 1927), Satz 172 


on p. 195. 
2R. Brauer, Uber die Darstellung der Drehungsgruppe durch Gruppen linearer 


Substitutionen (Dissertation, 1925), pp. 20-23. 
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ON THE MEAN VALUE OF r FOR KEPLERIAN SYSTEMS 
By Smmon PASTERNACK 
CALIFORNIA INSTITUTE OF TECHNOLOGY 


Communicated January 13, 1937 


I. Waller’ has proved, by means of the generating function for the 
Laguerre polynomial, that 


ee 
Se 





X the coefficient of (uv)*——! in 





2n (n+ 1)! 
(1 fee u)e+1(1 fs vet 
(1 — uy)2ts+8 : (1) 


where 1 is the total quantum number, / is the azimuthal quantum number 
and dp is the Bohr radius unit, a9 = h?/47°e?m. 
For s 2 —1, this yields the formula 


ee 
O \ag) Snatcy 4 (n—-1—1—k)! (n+D! 2) 
For s = —q,2 Sq S 21 + 2, we have 


n-l—1 
mR E(t) 
nay} 2n k=0 k 
(n-—1—1)!(n+1—q+1- )! 3 
(n-1L—1—2&)!(n+D! ss 
The latter expression is not suitable for computation in general, as the 
summation is from 0 to n — 1 — 1. 
J. H. Van Vleck? has derived an expression for r~, using Dirac’s ‘‘g- 


numbers.’’ A simpler expression may, however, be derived from (2) and 
(3) by use of the identity 


(21 —q)! 
(2/+q¢+ 1)! 


where 0 S gq S 21. 
Expression (2) may be put into the form 
a= (1) Ot a 
=\2z/ 2n (n+D! 

@ {(-o(-a +1)...(-9 +2 = Dn $141)... FF) 
feo ORL RI (—-n —1—g)(—n—1—q4+1)...(—n —l1—q+k-—1) 
nao\*-!1 (n +1+ q)! —-q¢ -g —n+1+1 
| 2 ee 2 » (5) 

2Z/ 2n (n+/)! 1, Seale 

















22)" it (4) 
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rt—2/A-l = ( 
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a, b,c 

e, f 

series with unit argument. Similarly expression (3) may be put into the 

form’ 

ate cme (et+i—q—1)!p.. Agtlati—s ie 

rw \oz/ on FDL te LL, ne - 2 +a +148 
(6) 


where AL | = F(a, b, c; e, f; 1) is the generalized hypergeometric 








To prove the identity (4), we use the following identity : 


gs Cee —b- orf) ae ee (7) 





ef | Me+f—b—-—or(f—a) ee +f —b—c 


Puttinga = —n+/1+1b=c= -ge=1,f = —n—I1— q+ 6, and 
then letting 6 —> 0, we get 


eed ee 
l, -n-l-q J (#414 Q)\(2l - 9)! 


lim ,;J@tlaqti,-m+/+1), 
— r en eee ir (8) 





The identity (4) follows at once from (5), (6) and (8). 
The computation of these mean values may be further simplified by 
using the following expressions, in which g 2 0: 


re ag (21 + 2g + 1)! ee en +n, alias 

















me =N\o7) on (+1)! 10+ 1,1 + */2 
r & wt 1 Qh + 2+ 1)! § 5, a + */s) [nt — C+ 1) 
. = 2n (21 + 1)! \ LI + 1). + */2) 
4 @= Vag + */2)(q + */2) [n? — E+ 1)"Mln? — C+ 2) ot 
21212 + 1) + */2)(0 + 2) + */2) 
(9) 
— nao\ * (21 + 2q + 2)! ee ee) 1+n,1+1 a 
== (27) (21 + 2)! PR pli ce 








= (=)" (21 + 2g + ae 4 4a + */)[n* — + 1?) 
Dae. (21 + 2)! 1!1!(2 + %/2)(1 + 2) 
(g — Valg + */2)(q + 8/2) [n? — (+ 17) fn? — C4+2)7) + 
2121(1 + */2)(l + 2)(1 + 5/2)(l + 3) 





+ 
(10) 


The expressions (9) and (10) follow from some properties of the poly- 
nomials® 
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l+2+m 
FM(z) = “| alii, saci “csi } (11) 


lrm+1 


It can be shown by means of the operational equation for these polynomials 
that 


m+t+1l+2m+i1-2 





—n, Nn + 1/5, 











2 : 2 
Fy,(s) = F (12) 
1, m+ i m+ 2 
2 2 
and 
1 a 
re —njn+y, @Pitam sy) | 

Fonti(Z) = “aut F cde een (18) 





wee | 


There is an identity® 


Ho b, 7 Pergrd — ard +6 — fri + bd — e) 


ef } Te-orf—ord+b—ald+b+c—fr(l +c—e) 


be—-a,l+bd-f ; 
ease eter (14) 


Puttinga = —g,b = —g,¢= —n+1+1e=1-—6,f=—-n-—l—gq+t 
6, and then letting 6 —> 0, we get 


a;* —q, —n +14 ‘ _ (= +q+ 1+)! 
1l,—-n—Il—q (+ 1)\n +14 Q)! 


ei year+ tt 
as } os 


A comparison of equation (15) with equations (5) and (11) indicates that 


— _ soi (—1)¢ (21 4 qd -f- 1)! Fa a 1(2n), 
2Z) “mm (Ql+)) 











(16) 


Equations (9) and (10) then follow immediately from equations (12) and 
(13). 
It can also be proved, by means of these generalized hypergeometric 
series, that 7 may be written in the form 
om @} 
ii (>) 1 feo + alll +1) + ad — DUC + 1) +2) + 
esl — 2) — IIL + 1) + 2)0 + 3) +...... (17) 
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where Cp, C1, C2, ...... are polynomials in m, and are odd or even functions 
of according as g is even or odd. 
There is a recurrence relation 


(n + 1)(n + m + 1)Feyi(2) = —Qn + leFe (2) + n(n — m)Fy-1(2). 




















(18) 
It follows then from equations (16) and (4) that, if g 2 1, 
2 x 242 
5 2q + 1 nao (2h + 1+ gl + 1 0S x (19) 
gq+1 Z qt+l 4Z? 
and 
1 (2 + ') Z 
a = _ r~a—2 
(+ 1+ 9/2) — g/2) \\q+1/ a 
a : 
i qg + 1 n’a)? mat (20) 


Formulas (4), (9) and (10) are the most useful for computational purposes. 


117. Waller, Zeit. Physik, 38, 644 (1926). 

2 J. H. Van Vleck, Proc. Roy. Soc. A, 143, 679 (1934). There is a slight error in his 
expression for r-* for 1=2, the coefficients being all 10 times too large. 

3 may be equated to zero except in the case g = 2). 

4W.N. Bailey, Generalized Hypergeometric Series (Cambridge Tract No. 32), p. 18. 

5 The polynomial F2(z) is discussed by H. Bateman, Tékohu Math. J., 37, 23 (1938). 
The general case will be treated by the author in a paper to be published soon. 

6 W.N. Bailey, Generalized Hypergeometric Series, (Cambridge Tract No. 32), p. 22. 


ON THE EQUATION OF DIFFUSION 
By Paut S. EPSTEIN 
NoRMAN BRIDGE LABORATORY OF PHysIcs, CALIFORNIA INSTITUTE OF TECHNOLOGY 


Communicated December 30, 1936 


1. Under the title ‘‘The Analytical Theory of Probability,” A. Kolmo- 
goroff! published a few years ago a mathematically elegant account of 
what might be called the general theory of diffusion. As the main result 
of his analysis, he sets up two differential equations for transition proba- 
bilities (see below). The first of them contains as its independent varia- 
bles the initial parameters of the transition probability and has an ana- 
lytical form resembling the ordinary diffusion equation. Unlike it, the 
second differential equation (depending on the final parameters) has a 
different form unfamiliar to the physicist. 
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The puzzling asymmetry of the two equations, as to whose physical 
reasons the author gives no hint, was, in our opinion, the main cause 
why the work of Kolmogoroff did not receive more attention. Unless 
these reasons are fully understood, an intelligent application of the theory 
to practical problems is very difficult. We think it, therefore, worth 
while to present in the following lines a discussion of the mutual connection 
between the two equations of Kolmogoroff and of their relation to other 
forms of the diffusion equation used in physics. 

2. The concept of the transition probability applies to systems com- 
posed of many identical elements. If one of these elements was at the 
time s in the state x, the probability of its being found at a later time / 
in the state y is called transition probability and is denoted by v(x, s; y, #). 
The states x and y are in general defined by several geometrical and 
dynamical coérdinates. Following Kolmogoroff, however, we shall 
restrict ourselves to the case of states depending on a single parameter, 
since the generalization for many parameters is obvious. Moreover, we 
shall take this single codérdinate as continuous and deal with intervals 
between x and x + Ax, etc. 

Kolmogoroff considers only the case of independent probabilities when 
v(s, x; y, t) is not influenced by the presence or absence of other particles 
in the same or neighboring states. Under this restriction this function 
must satisfy the two fundamental requirements of the theory of probability 


fo s; y,Hdy = 1, (1) 


foe s; y, t!)v(y, #; 2, dy = v(x, s; 2, t), (2) 


the integrals being extended over the whole region of variability of the 
parameter y. 

By mere mathematical transformations, these conditions lead directly 
to the equations of diffusion in Kolmogoroff’s form 


__ Ov(x, $5 yt) _ ) ; 0? : 
——— a(x, s) = v(x, s; y, t) + d(x, s) a v(x,s; y,8), (3) 

Ov(x,s; y,#) 9 F o? 2 
aie cane [a(y, t)v(x, s;y, #)] + or [d(y, #)o(x,s; y,#)]. (4) 


The meaning of the coefficients is the mean transition velocity of an 
element at the time ¢ and in the state y 


a(y, t) = limyeo - J (n — y)o(y, t; »,¢ + A)dn, (5) 
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and the mean spread of this velocity 
1 
b(y, t) = lima-o 5) f (n — y)*(y, t; t+ Addn. (6) 


Another type of probability which also has considerable physical interest 
is the distribution probability or density of distribution D(s, t; y) which 
is defined as the chance of finding an element in the state y at the time /, 
no matter where it came from. It is clear from the definition of the 
function v(x, s; y, t) that the density of distribution is connected with it 
in the following way 


D(s, t; y) = [06 $s; x)v(x, s; y, t)dx, (7) 


provided D(s, s; x) is the known density at the time s. Hence D(s, t; y) 
must also satisfy the equation (4). 


oD fe) 0? 
= “= [a(y, )D] + Dy [b(y, t)D]. (8) 

3. The most important case is the one in which the coefficients a(y) 
and b(y) are independent of time. The system is then called ‘‘homogeneous 
in time’’ and it is apparent from (5) and (6) that the transition probability 
can be represented as v(x, y, t-s), depending only on the elapsed time ¢-s. 
It will be sufficient for our purpose to consider this simpler case. In 
fact when we meet in the applications probabilities that are not homo- 
geneous in time they usually are also not independent and, therefore, 
outside the pale of Kolmogoroff’s theory. As the equation (8) shows, 
the main physical characteristic of these systems is that in them the den- 
sity of distribution D tends toward a state of equilibrium with the definite 
value Do(y). Supposing that the system has reached its equilibrium 
state, we can apply thermodynamical reasoning: For any interval of 
time the number of systems which move from the state x into the state y 
must be equal to the number of those which move in the opposite direction. 
This leads to the formula 


Do(x)u(x, y, t-s) = Doly)v(y, x, t-s), (9) 


as a direct expression of the second law of thermodynamics. 
For the function v(y, x, t-s), the first equation of diffusion (3) takes the 
simplified form 


ov(y, x, t-s) 


2) — aly) > 00,2, 1-8) + BO) 0,254). (10) 


Substituting from (9) v(y, x, t-s) = Do(x)v(x, y, t-s)/Do(y), we find 
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0 v(x, y, ts) _ fo E y, =) | »( El E y, =| 
ot Do(y) 4 oy Diy) a Oy? Dy(y) 


As a distribution density, Do(y) must be an integral (independent of 
time) of the same equation which is satisfied by v(x, y, t-s). We shall see 
that this requirement is fulfilled if we put 


C a 
Dy) = 5 exp Py dy, (11) 


where C isa constant. In fact, the equation for the transition probability 
v(x, y, t-s) becomes then 


ov re) 0? 
—=-— — [b ' 12 
ay = py Onl + 5, Ove] (12) 
while that for the equilibrium density can be written in the form 
re) 2 
0 = 2 [ay)Do) + = [60)Dah (13) 
oy oy" 


These equations are of the second Kolmogoroff type which is thus 
obtained from the first. The dsymmetry between the two types is, there- 
fore, shown to be an immediate consequence of the second law of thermo- 
dynamics. 

4. It is important to understand the physical nature of the coefficients 
a and 6 of the diffusion equation. According to the equation (5), a(y) is 
the mean transition velocity of all elements starting from the state y. 
At first sight, one should be inclined to think that systems in which this 
mean velocity vanishes must tend toward an equilibrium with uniform 
density of distribution. This is, however, not the case: Applying the 
second law of thermodynamics we found, in the preceding section,” that 
the equilibrium distribution Do(y) is not given by the general integral of 
the equation (13) but by the special integral (11). The condition for 


Dy being independent of y is, therefore, 2 exp( f : dy vo) oy = Oor 


db 
a= dy = 0. In other words, the mean transition velocity a contains 


a part (equal to db/dy) which has no tendency to change an existing uni- 
formity of distribution but the remainder a — db/dy is operative in dis- 
turbing it. This fact is in complete harmony with the considerations by 
which the equation of diffusion (8) is usually derived in physics. The 
argument starts from calculating the number of elements which move 
across a plane, y = const, per unit area and unit time. There are two 
classes of causes for this motion to be considered. Such causes as external 
forces or. temperature gradients belong to the first class and produce an 
effect of the form c,D simply proportional to the density. (For instance, 
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in the case of an external force F the coefficient c. = MF, where M is 
called the mobility). The causes of the second type are called forces of 
diffusion and produce an effect proportional to the negative gradient of 
the density —c®0D/dy. The accumulation of the elements in any point 
is, therefore, determined by the (negative) divergence of the total effect. 
This leads to the equation*® 


oD fe) re) oD 
— = — — (gD + 2(c oP) 14 
dt oy _ oy : oy (14) 
Comparing this with the form (8) we find 
b=, a — db/dy = q. (15) 


The expression for c is in agreement with the conclusions arrived at 
above. It shows that the form (14) of the equation of diffusion is a good 
guide to the physicist as it brings in evidence the quantities significant 
from his point of view. On the other hand, Kolmogoroff’s equations (3) 
and (4) are more general since they include the case of transition proba- 
bilities which are inhomogeneous in time. 

1A. Kolmogoroff, Mathematische Annalen, 104, 415 (1931). 

2 The substitution of the general integral of equation (13) for Do in (9) and (10) 
does not lead to the equation (12). Kolmogoroff claims the same results by imposing 
upon D, the condition lim Dp = 0 for y = +. That this argument is unconvincing, 
appears from the fact that it would exclude the uniform density. It is inconvenient 


to restrict oneself to finite systems. 
2 E.g., P. S. Epstein, Gerlands Beitr. Geophysik, 35, 154 (1932). 


ON THE STATISTICAL THEORY OF TURBULENCE 
By Tu. von KARMAN 
GUGGENHEIM AERONAUTICS LABORATORY, CALIFORNIA INSTITUTE OF TECHNOLOGY 


Communicated January 2, 1937 


G. I. Taylor! gave an important impetus to the statistical theory of 
turbulence by introducing the concept of “‘isotropic’’ turbulence, defined 
by the feature that the mean squares and mean products of the velocity 
components and of their derivatives are invariant with respect to rotation 
and reflection of the coérdinate axes. Taylor found that under the as- 
sumption of the isotropy, the squares and double products of the first de- 
rivatives of the velocity components can be expressed by one single correla- 


tion function R(y). This function is defined by the ratio “ where #, and 
U2 
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Ug are the simultaneous values of the x-components of the velocity fluctua- 
tion at two points at a distance y perpendicular to the x-axis. Taylor de- 
fines a characteristic length \ by writing 


: 2 
R(y) =1- = + higher terms in y. (1) 
2 
Obviously 2 ia ag (=) and Taylor shows that 
_BNdy#/, 9 

2 “2 
(*) = 22. (2) 

oy ty 


Taylor shows furthermore that the rate of decrease of the kinetic energy of 
2 
the unit volume by dissipation is equal to 15S (u = viscosity coefficient). 


He applies this result to the calculation of the decrease of the mean energy 
of the turbulent fluctuations downstream from a certain turbulence pro- 
ducing device such as a gridora mesh. Denoting the velocity of the main 
stream by U, he obtains the equation 


d-— 
i nal = ~10r = (3) 


where 4? and \? are both unknown functions of the distance x from the 
grid; vis the coefficient of kinematic viscosity. Taylor introduces a some- 
what arbitrary hypothesis connecting 4? and \? with the linear dimension 


1 
of the turbulence producing device, and tries to show that the reciprocal wa 


of the root mean square of the velocity fluctuation is always propor- 
tional to x. 

In this paper a second equation is established connecting 4? and \3, and 
it is shown that the decay of the turbulence is fully determined without 
further assumption by the correlation function R(y). Taylor’s result ap- 
pears as a special case corresponding to a certain type of the correlation 
function. 

We write the hydrodynamic equations in the form: 


~ + qVw — wVg = ww. (4) 


In this equation q is the velocity vector and w the vorticity vector. Mul- 
tiplying the three-component equations by é, 7, ¢ (where ¢, 7, ¢ are the 
components of w), respectively, and adding, the equation 
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D 
Di w? — w(w.V)¢ = w.V %w (5) 


bol = 


is readily obtained. 

We consider a fluid without mean velocity enclosed in a very large vessel 
and assume uniform, isotropic turbulence over the whole space considered. 
It is assumed that the vessel is so large that the influence of the walls can 
be neglected. Then taking mean values over time and the space consid- 


2 a a) 
ered, se will be replaced by ‘ w? and w.V*w by (V X w)?. The second 


term needs special consideration. Writing out the expression in codérdi- 
nates 


Ou Ov ow Ou , Ov 
(o.V)qg=8— +7 — + 2 — — a 
w.(w.V )g ee at Ot en (M4) 4 
ov , ow ow , Ou 
~~ ee — + — }.(6 
(= + oe) + te (* = 9) (6) 


It will be seen that the expression is built up as the sum of the products of 
the components of the tensor (ww) and of the components of the deforma- 
tion tensor. Now it can be shown that because of the isotropic feature of 
the correlations between the velocity components and their derivatives, 
expression (6) vanishes in mean value. The proof will be given elsewhere. 
As a matter of fact, any value of the expression (6) different from zero 
would mean that the vortex filaments had a permanent tendency to be 
stretched or compressed in the direction of the vorticity axis. With the 
expression (6) vanishing, equation (5) leads to the simple result: 


3 de? (2 <3] 
ca & = ioe oe 
2 dt : Oy ox @) 


and two similar equations obtained from (7) by cyclic permutation. De- 
fining a length X, by the relation 


2 ¢2 
() = 2% 
oy Ne 
which is analogous to (2), and computing the relations between squares and 
products of the derivatives of &, n, ¢, we obtain 
aire ¢2 
— = —10r-.. 8 
ne "x (8) 
Because of the isotropic feature of the system considered, £ can be ex- 


pressed by yu? and A?; namely £2 = 5 = Hence, equation (8) together with 
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Taylor’s equation (1) determines 4? as a function of the time, provided a 
relation between ) and X,, is known. 


We substitute in equation (3) . for U = and write 
dt dx 
u? 


— 
a w= — 1075; (9) 


Then we multiply (8) by u? and (9) by £ and subtract. Taking into ac- 
2 
count that A? = 5 =, we readily obtain 


d r? 
—r? = 10 te ~ 1), 10 
dt 2 ae 


. : ME 
Let us assume for the time being a constant value of the ratio. The 


physical significance of this assumption is discussed later. Then it follows 


cs ‘ ; ; , ? 
that \? is increasing at a constant rate with the time. Denoting (* ~ 1) 
by 6 








h? = do + 10frt (11) 
where Xj is the value of \2at¢ = 0. Then from (9) we obtain 
ma ui; 
uw? = (1 4 106» - (12) 
Xo 


denoting by x? the initial value of 4? at t = 0. : 
Applying equation (12) to Taylor’s problem, i.e., replacing ¢ by x/U, we 





obtain ; 
1 1 3 
SS = (1+ —]) . 13 
aa* dat oe si 
The symbols Vu, / uz denote root mean squares. Taylor’s result is 
: 1 x 
—= = —= + const. —. (14) 
Ve Vu U 


Our equation (13) gives this result as the first term of a series develop- 
ment. Obviously 
1 1 5B; x : ‘ 
— = —= + ——— — + higher terms in x. (15) 
Vu uw MVa U 
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If the linear law is correct for arbitrary large'values of x, 8 must be equal 
to 1/2. Now 6 can be expressed by the coefficients of the development of 
the correlation function R(y). We write 


2 
Ry) =1 “= eas ae 
First it can be shown that (=) = a(?) = 24by%. Then going 
oy? dy* /y=0 
through an analysis of the correlation between the various products of the 
second derivatives of u, v, w, similar to that performed by Taylor for 
the first derivatives,? we obtain 


By (tt 
~ 2 oy? 28 dy* i 


x 4+ lesa] ~~ 45 eal 
yo dy’? /y=0 
2 


: s r ° 

In order to estimate possible values of — the error function was first as- 
i _y 

sumed as correlation function, substituting R(y) = e ™. This leads to 


2 
‘ = 28/15 and 8 = 13/15 = 0.87. Taylor gives an example of the corre- 


w 
2 


lation function (see loc. cit., p. 45). Estimating b from the diagram >, = 


ow 


a as a function of x with values of B > 
Vue 
1/2, we obtain curves similar to those observed by H. L. Dryden. Unfor- 
tunately, no exact comparison between our theory and Dr. Dryden’s ex- 
perimental results is possible because apparently the decay of the inherent 
turbulence of the wind tunnel and the decay of the turbulence produced by 
the screens interfere with each other, and this interference modifies the 
law of decay, especially in the region of small intensity of turbulence in 
which the deviations from the linear law are significant. 

Taylor found that the constant factor in the second term at the right side 
of (14) is proportional to 1/M, where M is the linear dimension of the mesh 


which produces 1 the turbulence. Obviously, we obtain the same result by 


NV 72 


1.63 and 6 = 0.63. Plotting 





assuming proportional to M. Taylor made the arbitrary assump- 
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, » : : 
tion that “' remains constant and proportional to M along the whole 
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stream. In this new theory, it is sufficient to connect the initial values \; 
and 4? with the characteristic length of the turbulence producing device, 
the values of \? and 4? downstream are determined by differential equa- 
tions. It is felt that this represents a progress from the theoretical point of 
view. 

Considering the problem from a general point of view, it seems to the 
author that in an isotropic field the mean value of the turbulent fluctuations 
and the correlation function R(y) being given at a certain instant ¢ = 0, 
both y?2 and R are determined for all positive values of ¢. In other words 
we must be able to obtain a partial differential equation for R(y) with ¢ and 
y as independent variables. As an illustration of such a general theory, 
let us consider the case which can be described as ‘‘one-dimensional turbu- 
lence.’ Let us assume that the fluid has only one degree of freedom corre- 
sponding to a motion in the x-direction. We can imagine, for instance, that 
sheets of fluid extending parallel to the x-z plane are oscillating irregularly 
in the x-direction. The velocity is a function of y and ¢; however, it will be 
assumed that the oscillations of the fluid sheets are so fast and the ‘‘wave 
length’”’ of the oscillations so small that the mean value of 4? can be con- 
sidered as a slowly varying function of y and ¢. We consider the particular 
case in which x? is independent of y (homogeneous turbulence) and the 
correlation between the instantaneous velocities u, and “2 of two sheets at 
y = y, and y = ye is a function of the time ¢ and the relative distance 
n = Yo — ¥1. Let us now calculate the derivative of the product with re- 
spect to the time. Obviously 





ee Our OuUe 
— SS ane —. 16 
oy ama) —* + uy a (16) 
Now the hydrodynamic equations applied to the motion considered 
lead to 








2. 
Ou _ om (17) 
ot oy? 
and consequently 
6) — (= O72 ) 
— (uu) = v| — te + — 1 }. (18) 
ot ( 1U2 dy? 2 dye 1 
If the correlation function is denoted by R(n, ¢) 
um, = wR 
and 
Fy, Oe OR 


= . = 


U2 Uu : 
oy; os On? 
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Hence equation (18) can be written in the form 
2 


ee, — O'R 
7 (Ru?) = 2vy? —. (19) 


On? 
Applying this equation when 7 = 0 and taking into account that R(0) = 
1 and =) = 0, we obtain 


7=0 


d -— on or) 
= ‘(gy8) = Dp? | 20 
7 (u?) vu (= an (20) 
and eliminating x2 from (19) and (20) 
2 2 
OR, fOR_ (VR) | a) 
Ot On? On? 7=0 


The partial differential equation (21) is the fundamental equation of the 
problem considered. If R is given for ¢ = 0 as a function of 7, equation 
(21) determines the whole further development of the correlation and the 
mean value of the velocity fluctuation for all time. 

The most interesting case is that in which the shape of the correlation 
function does not change with time. If the shape of R remains similar, 


R(n, t) must be a function of the one dimensionless variable a, One of 
Ey. i 
these solutions is, for example, e 8. However, it is evident that in 
order to obtain such solutions, the initial shape itself must satisfy a cer- 
tain total differential equation. In general the shape of the correlation 
curve will vary with the time. The practical consequence of this result 
is the conclusion that simple rules like Taylor’s linear law for the decay 
of turbulence represent only approximations or special cases. 
Using for R(m) the development 


2 
Rn) = 1 — 5 + dn! (22) 
(=) ee 
On? pai 2 
and equation (20) reads 
5 ha te (23) 
dt ? 


which is analogous to our former equation (8) for the three-dimensional 
case. Substituting (22) in (21) and comparing the terms in 7?, we 
easily obtain 





Fo = YY 
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2 
o = 4y(6bA4 — 1) (24) 


which corresponds (putting 6b\? = *) to (10). From (23) and (24) it fol- 


‘wo 





lows that 
42 ue 
£(#') =— 4y (25) 
dt\\? NA 
. i ou — u? 
and introducing the vorticity fluctuation # = — Sy" w? = 2 rv 
y 
w? w? 
— = —4-—. 26 
at 2 (26) 


Hence, the equation (23) for the dissipation of energy and the equation 
(26) for the dissipation of vorticity appear as applications of the general 
partial differential equation (21), to small values of 7, i.e., to the immediate 
neighborhood of the point considered. In the special case that the corre- 
lation function preserves its shape throughout the whole motion, i.e., R is 
a function of r only, the ratio be remains constant; this was the assump- 

Vi w 
tion made earlier in this paper. Then the whole process of the decay of 
turbulence is determined by the initial values of u2, \ and X,, and we 
obtain equations analogous to (11) and (12). 

The equation corresponding to (21) for three-dimensional turbulence is 
given elsewhere. 

The author is indebted to C. B. Millikan and L. Howarth for discussions 
on the subject and for their assistance in calculations. 

1 Proc. Roy. Soc. A, 151, 421 (1935). 

2 Asa matter of fact, the analysis was carried out by means of a more general and more 


direct method than that used by Taylor. The new method which uses the conception 
of the ‘‘correlation tensor” will be published elsewhere. 
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INTRINSIC PROPERTIES OF LIGHT AND CORPUSCLES FROM 
DISTANT SOURCES 


By F. Zwicky 
NORMAN BRIDGE LABORATORY OF Puysics, CALIFORNIA INSTITUTE OF TECHNOLOGY 


Communicated December 31, 1936 


Spectra of distant nebulae, compared with those of neighboring nebulae, 
are shifted toward the red by amounts which increase with the distance 
which the light has traveled. This fact indicates that we are confronted 
with a phenomenon which involves history on a large scale, history either 
of the universe as a whole, or history of some of its individual parts. 

Scientifically speaking, history means the change in time of dimensionless 
ratios of significant physical quantities. Thus, on the relativistic inter- 
pretation of the redshift of light from nebulae the dimensionless ratio D/d 
between two lengths changes (increases) in time. Appropriate choices 
for these two lengths are Bohr’s characteristic length d = h?/41?me? as a 
supposedly fixed terrestrial measuring stick, and D = V/* where V is the 
volume which on the average contains one extragalactic nebula. However, 
many other interpretations of the redshift are possible. The assumption 
that history must be operative clearly suggests the necessity of an investi- 
gation of all dimensionless ratios between significant physical quantities. 
Only after this investigation has been completed will a final understanding 
of the redshift and other cosmic phenomena be possible. 

The following discussion will be concerned with the behavior of the 
most trivial dimensionless ratios only. A more general program may be 
outlined, but so far essential data for its realization are lacking. 


A. Nebular Light and Terrestrial Light—We shall designate one and 
the same property for nebular light and terrestrial light with py, and p,, 
respectively. The first step in our program consists in procuring data 
which will enable us to decide whether or not for a given property p the 
dimensionless ratio p,/p, is equal to unity or not. 

For instance, we may put py = Cy and py, = Cy where cy and Cp are 
the velocities of nebular light and terrestrial light, respectively, as measured 
simultaneously at the time of their arrival. Since no direct measurements 
are available we are forced to accept the values of cy obtained from the 
astronomical aberration as a legitimate substitute. Also, c, in the actual 
case refers to light from nearby stars instead of terrestrial light directly. 

The aberration of light from nebulae relative to the aberration of light 
from stars was measured by Strémberg and van Biesbroeck. Strémberg! 
has determined the positions of nebulae in the Ursa Major I cluster relative 
to the surrounding star field in the course of a year. The distance of the 
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Ursa Major I cluster is of the order of 72 X 10° light-years. Strémberg 
found that for light from these nebulae 


- Cy/Cp = 1 © 0.002. (1) 


This result indicates that if the redshift from nebulae is in any way 
related to a change of the velocity of light in time, this change can become 
apparent only if cy is measured relative to some standard other than the 
velocity of terrestrial light. This suggests that we investigate, as func- 
tions of time, ratios of the type c/v where v is any significant prototype of 
a velocity other than c,. For instance we might choose v = 2ze?/h, 
where v is the “‘velocity”’ of an electron in the first ‘‘orbit’”’ of a hydrogen 
atom. Such a dependence of c/v on time might exist if the physical 
properties of the universe as a whole were subject to history. Although a 
few experiments have been carried out with some bearing on the above 
problem, no direct data are available from which one might deduce the 
behavior of c/v over long intervals of time. We therefore drop this matter 
for the present. 

Our next step calls for a comparison of the interaction with matter of 
nebular light and terrestrial light. 

We define the wave-length \ of a monochromatic beam of light by the 
angle 6 which the diffracted beam forms with the original beam of normal 
incidence on a plane grating whose spacing is A. The definition equation 
is 

\= Asin 0. (2) 


Assuming that A/d is independent of time, wave-lengths, through the 
observation of diffraction patterns, are measured in terms of the funda- 
mental atomic length d = h?/4r*me?. The frequency vy is not measured 
but is introduced by definition as v = c/\. Absorption measurements 
and the photoelectric effect show that light of the wave-length \ acts in 
many ways as if it were composed of photons characterized by individual 
energies « = hv = hc/d, where Planck’s constant / is presumed to be a 
universal constant. If this is correct then we should have h,/h, = 1. 
Observations to test this point have been made by W. S. Adams and 
M. L. Humason.? These observers have determined the redshift relative 
to a terrestrial source of light from the extragalactic nebula N. G. C. 4151 
both with a grating and a prism spectrograph. The refraction measure- 
ments yielded an average apparent velocity of recession 1g = 958 km./sec., 
whereas the grating gave vg = 962 km./sec. Assuming that only an 
insignificant part of v is due to the peculiar motion of the nebula it follows 
from these observations that 


(he) y/ (he) = 1 * 0.02. (3) 
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In combination with (1) we obtain 
hy/hp = 1 © 0.02 (4) 


and we may conclude from (1) and (4) that nebular light and terrestrial 
light possess intrinsically identical: physical properties. 

In view of the fact that the nebula used in the above observations might 
accidentally have a peculiar motion with a large component of recession 
it will be advisable to repeat the observation with some more distant 
nebula and to use an absorption test (rare earth filter) instead of the re- 
fraction test. If this test gives the same result as the one made by Adams 
and Humason we will have ascertained that nebular light at the time of 
arrival has intrinsically the same properties as terrestrial light. Assuming 
this to be correct, the existence of a redshift simply means that certain 
characteristic features, such as maxima and minima of intensity, or emis- 
sion and absorption bands in the spectra of distant nebulae appear to be 
displaced toward the red as compared with what we think to be comparable 
features in the spectrum of a nearby nebula or in the spectra of appro- 
priately chosen stars and terrestrial sources. Continuing our scrutiny 
of apparently trivial facts, the question then arises as to what we mean by 
comparable features. The analysis of this question which leads to the 
investigation of a number of fundamental dimensionless ratios will be 
given in another paper. 

Before abandoning this discussion of the intrinsic properties of light I 
sketch briefly a related problem which may find its solution through 
astronomical observations. Certain effects have recently come to our 
knowledge which suggest that the differential equations governing the 
propagation of light through empty space are not strictly linear. Non- 
linearity of these equations results in a dependence of the velocity of light 
on frequency. One reason for the existence of slight deviations from the 
superposition principle of light lies in the potential possibility of the forma- 
tion of pairs of positive and negative electrons by interacting photons. 
This interaction necessitates a generalization of Maxwell’s field equations 
through the introduction of non-linear terms in the field strengths.’ As 
a result, light traveling through space which is free of matter but filled 
with radiation will have a velocity depending on its frequency. In addi- 
tion similar effects arise from the gravitational interaction of light with 
light and matter. Although the effects to be expected are in all probability 
small the possibility of an experimental test may be kept in mind. 

An obvious way of investigating effects of the kind mentioned lies in 
the observation of light signals of different frequency which have traveled 
through space for a long time. If we knew of any signals which have 
started simultaneously from a very distant source the dependence of the 
velocity on frequency could perhaps be demonstrated by checking up on 
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the times of arrival of these signals. A unique case to perform this test 
is afforded by the observation of distant nova outbursts. For such an 
outburst we may safely assume that photons in the various emission 
lines of, say, hydrogen have started out simultaneously on their long 
journey. We must therefore attempt to determine whether or not the 
intervals in which the different hydrogen lines flare up in the spectrum of 
a distant super-nova are zero or not. With the present telescopic equip- 
ment it is probably possible to find super-novae to distances as great as 
10° light-years. With reasonable luck it should be possible to detect 
differences in the time of arrival of various emission lines amounting to as 
little as a few days. Observations of this kind will therefore enable us in 
the most favorable cases to determine the ratio of the velocity of violet 
light to that of red light traveling through internebular space with a relative 
accuracy of 10~-". 

There already are available some data which might have a bearing on 
our problem. These data refer to a difficulty which has arisen with 
respect to the absolute magnitudes and the shape of the light curves of 
common novae in the Andromeda nebula (Messier 31). About one hundred 
uovae have been found in this galaxy, a number which would seem suffi- 
ciently large to constitute a fair collection and to make possible a statistical 
comparison with the characteristics of galactic novae. This comparison 
results in two serious discrepancies. 

In the first place no common nova in M31 has been observed whose 
absolute magnitude is brighter than M = —6.7. Galactic novae on the 
other hand often are as bright as M = —8 oreven M = —9. This dis- 
crepancy was pointed out to me by Dr. W. Baade of the Mt. Wilson 
Observatory. 

In the second place the light curves of the novae in M31 are on the 
average considerably flattened as compared with the light curves of similar 
galactic novae. In particular, the familiar flash or peak which forms the 
maximum of the light curves of many galactic novae has apparently never 
been exhibited by any of the novae in M31. 

Both of these discrepancies could readily be understood if the velocity 
of red light and violet light differed by about 10—* to 10-? km./sec. 

I have also tried to check up on the light curves of Cepheid variables in 
nebulae. Extragalactic variables of a given period should have flatter 
light curves than galactic variables of the same period if the velocity of 
light depends on frequency as strongly as conjectured above. The data 
available at the present unfortunately are not accurate enough to test 
this point. 

Another startling fact was observed by Sinclair Smith’ of the Mt. Wilson 
Observatory in the spectrum of the recent super-nova in N. G. C. 4273, a 
member nebula of the Virgo cluster. The spectrum of this nova, contrary 
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to all expectations, did not contain any emission lines in the ultra-violet 
region between 1 = 3100 A and A = 3600 A. In fact the intensity of 
light in this region which might be ascribed to the super-nova is indis- 
tinguishable from zero. This fact again can be explained on the assump- 
tion that the ultra-violet light from the super-nova mentioned did not 
arrive simultaneously with the visible light and that it hence was missed. 
Alternative explanations, however, are still possible. It will be of great 
interest to discover super-novae in as early stages as possible and to care- 
fully analyze their spectra. 

Provided that cosmic rays originate in novae* the study of fluctuations 
in the intensity of these rays will be of particular interest in regard to 
the problems formulated in this paper. So far, however, the available 
data are insufficient to draw any definite conclusions. 

B. Corpuscular Rays.—Common novae, super-novae, Wolf-Rayet 
stars, etc., eject into interstellar space energetic corpuscular rays of all 
descriptions. On the other hand, investigations of the polar and the 
non-polar aurora, as well as the cosmic ray phenomena, indicate that 
corpuscular rays are continually impinging on the earth’s atmosphere 
from extraterrestrial regions. To correlate these phenomena with the 
happenings on unstable stellar objects is very suggestive. It is also sug- 
gestive to extend the inquiry proposed in this paper from photons to 
corpuscular rays and to ask whether or not extraterrestrial corpuscles 
exhibit the same intrinsic properties as terrestrial corpuscles. To answer 
this question one must investigate dimensionless ratios of the type py/pz, 
where p = e, m, uy, etc., are the charge, mass, characteristic wave-length, 
magnetic moment of a corpuscle. Little evidence is available at the 
present which would enable us to determine any of the significant ratios 
Py/pr. It will be well, however, to keep in mind the eventuality of 
values p,/p,; # 1 for the interpretation of cosmic ray phenomena in 
order to prove or disprove the so far tacitly made assumption of the 
uniformity of all happenings in the observable parts of the universe. 

The program submitted in this paper was outlined several years ago 
and was discussed on various occasions with those astronomers at the 
Mt. Wilson Observatory who have made the investigation of distant 
nebulae their specialty. I am particularly indebted to Drs. W. Baade, 
E. Hubble and M. L. Humason for codperating in the attempt to test 
some of the points of this program. 

1G. Strémberg, Pub. Astr. Soc. Pac., 43, 266 (1981). 

2 W. S. Adams and M. L. Humason, Pub. Astr. Soc. Pac., 48, 107 (1936). 

* I am indebted to Dr. Sinclair Smith for private communication of his data. 


4 See W. Baade and F. Zwicky, these PROCEEDINGS, 20, 259 (1934). 
5 Ruler and Kockel, Naturwiss., 23, 246 (1935). 





5. VoL. 23, 1937 ASTRONOMY: D. HOFFLEIT 111 


ot SPECTROSCOPIC ABSOLUTE MAGNITUDES OF THREE 
of HUNDRED AND SEVENTY SOUTHERN STARS 


. By Dorrit HoFFieir 
>t HARVARD COLLEGE OBSERVATORY, CAMBRIDGE, MASSACHUSETTS 


Communicated December 29, 1936 


me The Mount Wilson Observatory has published spectroscopic absolute 

magnitudes for more than four thousand stars north of declination —26°;! 
- the Harvard Observatory, in obtaining comparable data for stars south of 
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FIGURE 1 
A composite Russell diagtam. Dots represent Harvard 
data; open circles, Mount Wilson data. Corresponding 
points are joined by straight lines. Ordinates are absolute 
magnitudes, abscissae are spectral classes. 


—20°, is extending the survey of spectroscopic parallaxes to the South 
Pole. The spectra were made with the 13-inch Boyden telescope and a 
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singe objective prism (dispersion = 22 mm. between Hf and He). Absolute 
magnitudes have now been determined for 370 stars with spectral types 
later than A5 and brighter than the fifth apparent visual magnitude. The 
present note summarizes the principal results; the detailed investigation 
is to be published later in the Harvard Annals. 

1. For the determination of absolute magnitudes an accurate classi- 
fication of the spectra is of importance. On account of the well-known gaps 
in the Draper system of classification it was necessary to sub-classify the 
spectra for the 370 stars. For this purpose a few selected spectra were 
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FIGURE 2 





Dependence of M;,-M; on spectral class. Dots represent trigonometric parallaxes 
20°020; opencircles, trigonometric parallaxes <0"020. Diagram at left gives Harvard 
data; at right, Mount Wilson data for bright stars in first 8 hours of Right Ascension 
and the stars common to the two lists. 


used as standards for comparison. The revised Harvard system differs 
slightly in scale from the system of classification used at Mount Wilson, 
viz., the Harvard class F5 corresponds to the Mount Wilson F4, and the 
Harvard K5 to the Mount Wilson K6 or MO. 

2. Absolute magnitudes were derived from the ratios in intensity of 
various pairs of lines. The number of pairs available within any one 
spectral class varied from four to ten. For calibrating the intensity ratios 
three different methods were used based on (a) spectroscopic absolute 
magnitudes determined at Mount Wilson for the 140 stars common to the 
Harvard and the Mount Wilson lists, (b) trigonometric parallaxes and (c) 
proper motions. The zero points of the mean calibration curves were 
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corrected to bring the final system of absolute magnitudes into close agree- 
ment with the Mount Wilson system. 

3. From the agreement between the magnitudes derived from the 
various pairs of lines, the average probable error of a final mean absolute 
magnitude is found to be 015. The average difference between the 
magnitudes determined from the two important pairs 4077:4071 and 
4216:4250 is 0°65. These two errors are of the same order as the cor- 
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FIGURE 3 
Dependence of M,)-M; on trigonometric parallax. Dots indi- 
cate high galactic stars; open circles, stars within 15° of galactic 
equator. The five discrepant points in the upper diagram all 
represent c-stars. 


responding errors of the northern stars studied at the Dominion Astro- 
physical Observatory ;? other available catalogs are not published in suffi- 
cient detail for similar comparison. 

4. The differences between the recently determined absolute magnitudes 
for southern stars and the absolute magnitudes for the same stars deter- 
mined elsewhere are summarized below: 
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Sources or M DIFFERENCE No. Stars 
AVERAGE SYSTEMATIC 
Mt. Wilson minus Hoffleit +()M40 +0M04 137 
Rimmer? minus Hoffleit +0 .66 +0.21 103 
H.C. 228, 232 minus Hoffleit +0 .82 —0.72 46 
H.C. 243, 246 minus Hoffleit +0 .50 —0.18 162 


The values in Harvard Circulars 228 and 232 are based on the system of 
Mount Wilson Contribution 199, those in Circulars 243 and 246 on proper 
motion data. 

5. A comparison of the Harvard and the Mount Wilson systems of 
absolute magnitudes is given in the composite Russell-diagram (see Fig. 1). 
Only the stars common to the lists of both observatories are included. 
The systematic difference in spectral class, already noted above, is the 
most conspicuous feature of the comparison. For the supergiants and 
the few intermediate giant stars the mean differences in absolute magnitude 
are large because of insufficient calibration material for such stars at 
Harvard. The mean differences in absolute magnitude between the 
Harvard and the Mount Wilson systems show no dependence on spectral 
class. 

6. The spectroscopic absolute magnitudes (M,,) have been compared 
with absolute magnitudes computed from trigonometric parallaxes (M,) 
taken from the Yale Catalogue (1935). The dispersion in the differences 
M,, — M, increases with advancing spectral type (Fig. 2). The average 
difference M,, — M, does not vary with spectral class if no trigonometric 
parallaxes, large or small, are rejected; but if we reject parallaxes with 
errors greater than half the parallaxes themselves (the open circles in 
Fig. 2), then the mean difference M,, — M, increases negatively with the 
later types. 

7. For the stars with small trigonometric parallaxes, the differences 
M,, — M, show a pronounced dependence on the size of the parallax 
(Fig. 3). The correlation may be attributed in part to systematic errors 
in the trigonometric parallaxes: for example, a constant error of —0"005 
in the trigonometric parallaxes would, in the absence of any other sources 
for systematic errors, account for the run in the observed differences 
Msp — M;. 

1 Mt. Wilson Contr., No. 511; Astrophys. J., 81, 57 (1985). 


2 Pub. D. A. O., 3, No. 1 (1924). 
3 Mem. R. A. S., 62, 113 (1923); 64, 1 (1924). 
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CHANGES OF PERIOD IN VARIABLE STARS OF LONG PERIOD 
By T. E. STERNE AND L. CAMPBELL 
Tue HARVARD COLLEGE OBSERVATORY 


Communicated December 23, 1936 


One of us has investigated elsewhere! the mathematical theory of the 
accuracy of periods determined from dates of maximum, or minimum, or 
mth magnitude increasing or diminishing, of a variable star. The existence 
of cumulative errors? in the periods of long-period variables was shown to 
provide a simple explanation for the major part of their observed irregu- 
larities of period, as revealed for instance by O-C diagrams,’ and plots 
of cycles against epoch number. The existence of cumulative errors was 
shown to be in itself sufficient to cause an erratic behavior in O—-C diagrams, 
so that such diagrams did not by themselves furnish safe evidence for 
changes of true, or mean, period. Further, it was shown that a value of 
the mean period, obtained from observations extending over any finite 
interval of time, would be liable to errors of sampling, so that the mere 
existence of a scatter among values of the mean period obtained from 
differing intervals of time was not in itself evidence for a changing true 
or mean period. 

Before the appearance of the previously mentioned investigations, it 
had been the common practice among workers in the field of variable stars 
to identify irregularities in O-C diagrams with “‘changes of period” in the 
stars. Thus, of the 566 variables with known long periods listed‘ in the 
Katalog der Elemente des Lichtwechsels fir 1687 veranderliche Sterne, sine 
terms were assigned to 85. Many of the reported “‘changes of period”’ 
were possibly of the nature of mere statistical fluctuations among cycles 
having a natural dispersion. It therefore became eminently desirable that 
the matter of changes of period should be reéxamined at length, in an 
impartial and systematic fashion. 

The authors have examined systematically, for changes of true (or mean) 
period, 377 variable stars with periods usually in excess of 100 days. The 
stars include nearly all of the best-observed long-period variables. Ap- 
parent “changes” of period, that can reasonably be ascribed entirely to 
statistical fluctuations, have been carefully distinguished from systematic 
changes of period, here called changes of true or mean period. The material 
examined was obtained from 


(1) observations by members of the American and French variable star 
associations, and 


(2) earlier observations, referred to in (a) the extensive card catalogue 
bibliography of variable stars, compiled by Miss A. J. Cannon at Harvard, 
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and in (b) the bibliography contained in the Geschichte und Literatur des 
Lichtwechsels, edited by the Astronomische Gesellschaft. 

The method of study employed the x? test of Karl Pearson. More 
complete details will be published in the Annals of the Harvard College 
Observatory. We have obtained the following results. 

1. The SS Cygni type stars, SS Cygni, U Geminorum and SS Aurigae, 
are found to have sinuously changing true periods. 

2. The Z Camelopardalis type stars, RX Andromedae and Z Camel- 
opardalis, show evidence of change in the true periods of the 120 in- 
creasing phases, but show no good evidence of change in the true periods 
of the 12"0 decreasing phases. The existence of sinuous changes in the 
widths of maxima would constitute a possible explanation of the difference 
between the results obtained from the two phases. 

3. Twenty of the 372 remaining variables are ‘‘suspected’’ of having 
changing true or mean periods, but of these only about five can be stated 
with any high degree of assurance, on the evidence, to have changing true 
periods. Although it is not certain which five these are, the most reason- 
able choices are R Hydrae, R Aquilae, R Cancri, U Bodtis and S Serpentis. 
The first two have probably had diminishing true periods, the remaining 
three sinuously changing true periods. 

4. The small percentage of the stars which can now be stated with 
assurance to have changing true periods arises probably in part from the 
insufficient duration and number of observations. However, among 52 of 
the stars, not of the SS Cygni or Z Camelopardalis types, observed since 
before J. D. 2,400,000 and therefore longest observed, only about four can 
be stated with assurance to have changing periods. These 52 stars have 
been observed, on the average, for over ninety-four years. 

5. Most of the previously reported ‘changes of period’ among long- 
period variables can be reasonably interpreted as mere statistical fluctua- 
tions. 

6. There is no good evidence for any secular increase or diminution of 
period, of an evolutionary nature, among the long-period variables as a 
whole. 

7. We have obtained an accurate list of periods. 

In closing, the authors wish clearly to point out that no conclusions 
as to the proportion of variables with changing periods, reached by this 
or any other investigation, can be regarded as final. An increase in the 
interval of time covered by the observations, and in their number, must 
render any tests more sensitive to continuing changes of period. All that 
we have been able to do has been to find out how many of the 377 variables 
which we have studied may be said with assurance, on the basis of observa- 
tional material now available, to have changed their periods in a manner 
which cannot be entirely attributed to mere statistical fluctuations. 
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1T. E. Sterne, Harvard Circular Nos. 386, 387 (1934); Popular Astronomy, 42, 558 
(1934). 

2 The existence of cumulative errors in long-period variables had previously been 
inferred by Eddington and Plakidis, Monthly Notices, 90,65 (1929). See also Plakidis, 
Ibid., 92, 460 (1982); 93, 373 (1933). 

3 O-C diagrams are plots, against epoch number, of the difference, observed date 
minus computed date, of the phase selected for study. 

* Geschichte und Literatur des Lichtwechsels, 3 (1922). 


THE SYNTHESIS OF GAMMA-METHYL-5,6-CYCLOPENTENO.- 
RETENE, A COMPOUND STRUCTURALLY RELATED TO “DIELS’ 
HYDROCARBON”! 


By Davip E. ADELSON? AND MARSTON TAYLOR BOGERT 
ORGANIC LABORATORIES, DEPT. OF CHEMISTRY, COLUMBIA UNIVERSITY 


Communicated January 12, 1937 


By fusion with selenium, the so-called ‘‘Diels’ Hydrocarbon’”’ (I) (gamma- 
methyl-1,2-cyclopentenophenanthrene) has been obtained from such 
interesting biologically and chemically important compounds as cho- 
lesterol, ergosterol, lumisterol, isopyrovitamin, pyrocalciferol, cholic acid, 
strophanthidin, gitogenin, anhydro-uzarigenin, sarsasapogenin and pseudo- 
bufotalin. 

Retene (II), likewise a phenanthrene derivative, is formed in the de- 
structive distillation of coniferous woods, or of rosin, due to the breakdown 
of the abietic acid present, and is found in large amount in the higher 
boiling fraction of the resultant tar. When abietic acid itself is heated 
with either sulfur or selenium, retene is the chief product. Further, 
retene carries the methyl and isopropyl groups, two radicals which are 
frequently found associated in natural products (cymene, terpenes, cam- 
phors, etc.), and which seem to have some significant function in nature. 
The attachment of the isopropyl group in the retene molecule, for example, 
exhibits the same carbon skeleton as is found in that important natural 


building unit isoprene, viz.: 
~~ Cc 
Cc 


It seemed to us, therefore, desirable to attempt the synthesis from 
retene of products akin to the Diels’ hydrocarbon and related compounds, 
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in the hope of leading ultimately to the establishment of closer chemical 
and biological connections between the resin acids and such other note- 
worthy phenanthrene derivatives as the sterols, bile acids, hormones, 
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(V) Carbon skeleton common (VI) 6-Acetylretene 


to (III) and (IV) 


cardiotonic drugs, toad poisons, morphine alkaloids and the like; and of 
throwing additional light upon the relation between chemical constitution 
and physiological action in these interesting compounds. 

This applies also to the cancerigenic and estrogenic properties already 
recorded for certain phenanthrene derivatives. Of those so far investi- 
gated from this standpoint, the most potent cancerigenic hydrocarbon is 
methylcholanthrene (IIT). 

The new hydrocarbon synthesized by us is a gamma-methyl-5,6-cyclo- 
pentenoretene (IV) and, as will be seen by a comparison of the structural 
formulas, differs from the Diels’ hydrocarbon (I) (which is neither canceri- 
genic nor estrogenic) in the location of the methylcyclopenteno ring and in 
the presence of the additional methyl and isopropyl groups. From the 
methylcholanthrene (III), it likewise varies in the position of the cyclo- 
penteno cycle, and in lacking the 4th hexacarbocyclic system. The carbon 
skeleton common to both (III) and (IV) is shown in (V). Some of our 
good friends in the medical profession are generously codperating by ex- 
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amining this new product pharmacologically, to determine whether or not 
it possesses cancerigenic properties. 


The steps included in the synthesis of (IV) from 6-acetylretene (VI) are 
shown below, in which CisH;7 signifies the retyl residue: 


CisHi7COCH; + BrCH,COOR + Zn — 
OH 


CisHiyC(CH3)CH2,COOR + KOH + MeOH —> 
OH 


CisHivC(CHs3)CH2COOH + (CH;CO),0 + CH;COONa —> 
CisHi7C(CHs): CHCOOH + Na-Hg —> 
CisHi7CH(CH;)CH2COOH + PCI, —> 
CisHi7CH(CH;)CH2COCI + AICI; —> 


CH(CHs) 
CiHis CH, + Zn-Hg + HC] —> 
CO 
CH(CHs) 
CisHis CH:; short white needles, m. p. 74.5-75.5° (corr.); soluble 
CHy 


in the common organic solvents and in vegetable oils. Absorption spectra, 
x-ray crystal measurements and other physical data, as well as other de- 
tails concerning the above products and certain of their derivatives, will 
appear in full elsewhere. 

Similar syntheses leading to benzoretenes are nearing completion and 
will also be published shortly. 

The intense activity in this field at present, both here and abroad, has 
made it seem wise to issue this preliminary communication without further 
delay. 


1 Reported to the Symposium on Cancer, Am. Assoc. Adv. Sci., Atlantic City Meeting, 
Dec. 30 (1936). This is a preliminary communication; full experimental details will 
appear in one of the chemical journals. 

2 NATIONAL RESEARCH COUNCIL FELLOW. 











120 CHEMISTRY: W. D. HARKINS Proc. N. A. S. 


THE INTERMEDIATE NUCLEUS IN THE DISINTEGRATIVE.- 
SYNTHESIS OF ATOMIC NUCLEI: DISINTEGRATION IN STEPS* 


By WILutraM D. Harxins** 


GEORGE HERBERT JONES CHEMICAL LABORATORY, UNIVERSITY OF CHICAGO, AND BAKER 
CHEMISTRY LABORATORY, CORNELL UNIVERSITY 


Communicated December 21, 1936 


1. Introduction: Nuclear Reactions——This paper deals with the most 
fundamental problem of nuclear reactions. While much has been written 
concerning the artificial disintegration of atoms, very little attention has 
been paid to their synthesis. For a considerable number of years the 
writer has used the point of view that the primary process is always the 
synthesis of a new nucleus, which disintegrates only to get rid of the excess 
of energy, over that of a stable state, which has been given to it by the 
synthesis. Thus the first nuclear reaction to be studied was that in which 
a fast helium atom, with a velocity greater than that which corresponds to 
an energy of more than three million electron volts, collides with the nu- 
cleus of a nitrogen atom ina gas. The reaction which occurs is the synthe- 
sis of an atom of fluorine: 

oni4 +. $He* = 0 F*16 ; (1) 
Here for nitrogen (VV) the 7 represents the atomic number or nuclear posi- 
tive charge, the 14 gives the atomic mass number (A) and the 0 represents 


the isotopic number J, 
I=A-—2Z (2) 


which is very important in the classification of atomic nuclei. 
The composition of any nucleus is represented by the formula of 
Harkins! (1922) adopted later by Heisenberg? which is 


(pn) gn, (3) 


in which is a proton, and m a neutron. 
The excited fluorine nucleus has a short life and may disintegrate by 


either of the two methods: 


SFris ae 3017 +4 1H (4) 
or 

o F*i8 per Fer 4. yn} (5) 
followed by 

“LRH _, JO" 4 —te0, (5’) 


2. Theories of Nuclear Disintegration.—There have been two theories of 
the mechanism of the nuclear disintegration of atomic nuclei, as follows: 
(1) Theory of the Intermediate Nucleus. 
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This type of mechanism, as suggested by the writer and outlined above, 
was opposed by nuclear theorists prior to February 1936. 

(2) Theory of Direct Disintegration, or Non-Existence of the Interme- 
diate Nucleus. 

This latter theory was supported by all quantum theorists interested in 
the nucleus. 

The idea that a nuclear reaction occurs in at least two steps (Theory 1) 
and not as a single step, was expressed by the title of an early paper* ‘“‘The 








‘2m 
146° 
128° 22mm 
| 86° 
| FIGURE 1 
! This figure gives the actual dimensions for the 
{oF mm capture of a neutron (dotted line) by a nitrogen 
| nucleus of mass 14 (upper end of dotted line) to 
| give a nitrogen atom of mass 15, which disinte- 
| grated into a boron 11 atom with a track 2.2 mm. 
| long and a helium atom (track 12 mm. long). 
46mm 
| The neutron had a velocity of 1.89 X 10% cm. 
Gass fs 
[ per second, and an energy of 3.00 million elec- 
| tron volts, and passed in a straight line through 
| 9 millimeters of lead, 16 mm. of glass and 64 mm. 
of air. The velocity corresponds to 14700 miles 

| 9mm per second. This was the first neutron found in 
| Lead America. 
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Synthesis and Disintegration of Atoms.’’ This expression was later 
shortened to “the disintegrative-synthesis of atoms.” 

An independent type of theory” was favored by nuclear quantum theorists. 
According to this point of view a disintegration is always direct, without 
an intermediate nucleus, but may occur in two ways: 


(a) By capture of the projectile, 
For example 
N'“ + Het — O" + . (6) 

According to this theory the intermediate nucleus has no life and no exis- 
tence: There is only a trading of parts. 

(b) Disintegration without capture: “‘pure disintegration.” 

That is, a nucleus, P, used as a projectile, strikes a nucleus 7, used as a 

target, and rebounds, while T splits in two (or possibly more) parts: 
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T—> C+D. (7) 


This point of view presented by Gamow, and supported by Chadwick and 
others, was that generally adopted, as corresponding to one type of 
disintegration. 

That neutrons, at least, can give ‘‘non-capture”’ or pure nuclear disin- 
tegrations seemed to be demonstrated by the work of Feather‘ at Cam- 
bridge, and may be expressed in his own words: 

“About 130 cases of interaction between a neutron and a nitrogen nu- 
cleus have been observed; of these about 30 resulted in disintegration, 
more than half of the latter without capture of the neutron. This is very 
different from the results obtained under a-particle bombardment, where 
elastic collisions, resulting in measurable spurs in an expansion chamber, 
outnumber inelastic (disintegration) collisions by a factor of the order of 
1000:1. Moreover, although the possibility of non-capture disintegration 

by a-particles has frequently 


D D been considered, unexception- 
. 7? . able evidence for its occurrence 
c “ é 5 Pa ———K has yet to be obtained. The 
A 
A 


former of these points of dif- 
ference is certainly to be ascribed 
A to the different extent of the 
external fields of the two par- 








/ 
a - ticles, that of the neutron fall- 
db ef } ing off very rapidly to become 
FIGURE 2 © already inappreciable at a few 

Capture. Non-capture. diameters distance; it is quite 


possible that further investiga- 
tion will exhibit the latter difference, also, as a direct result of the same 
circumstances.”’ 

The disintegration “‘by capture’ is illustrated by figure 1, in which the 
dotted line indicates the invisible track of the neutron. Figure 2c illus- 
trates the disintegration by “‘non-capture”’ as assumed by Feather and by 
Chadwick. Thus the reaction must occur by non-capture if the neutron 
comes directly from the source. It was assumed, however, by Harkins, 
Gans and Newson,* that an intermediate F'* is formed, by the action of a 
scattered neutron (Fig. 2b). According to them the reaction is 


74 + ont — 7N 15, (8) 
1N15* _, (Bu + dHet. (8’) 
3. Rapid Shift of Theory.—In.a paper written early in 1933, the writer,§ 


who believed that the synthesis of an intermediate nucleus is an essential 
preliminary to a nuclear disintegration, assumed that all of these disinte- 
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grations which seem to be obviously by non-capture, occur actually by cap- 
ture, but that in each case the neutron has been deflected (scattered) by an- 
other nucleus. (Fig. 2b). 

It seemed that an essential step in the advancement of the theory of the 
formation of the intermediate nucleus, was to secure evidence against the 
non-capture disintegrations. To obtain this Harkins and Gans,® in 1933- 
34 developed the (relativistic) mechanics for a non-capture disintegration. 
This had not been done previously, and it was at once obvious that all of 
the apparently non-capture disintegrations which had been obtained were 
actually by capture. Thus the velocities for the projectile given by these 
equations are much higher in general than those necessary if by capture, 
as is apparent from the equations as given below. The velocities were found 
to be in general much higher than any possible velocities for the actual 
projectiles, as may be deduced from the equations given below. The rela- 
tivity equation for the velocity of the projectile is: 


(Rome + Rpmpvp)? — c(keme + kRpmp — mz)? + 
— 2c?k 4m (Rome + kpmp — mz) 
Z 2k ,m4(Romgtg + kpmypvp) cos a 





while the similar equation for not too high velocities is 


- M? + 2m,E, + Ep + (me + mp — mz) 
2m,M cos a 





Va 


in which M represents momentum, E is kinetic energy, m mass, v velocity, 
c the velocity of light and A, B, C and D are the nuclei involved. If a 
y-ray is emitted the expression between brackets in equation (16) contains 
in addition a term Ey. 

4. Theory of the Intermediate Nucleus. (Fig. 3).—As an example of the 
use of the idea of the intermediate nucleus, Harkins, Gans and Newson® 
consider that in the production of neutrons by bombardment of atomic 
nuclei by a-particle (}He*) the first step is an increase in the atomic number 
to a value 2 higher than that of the target atom, and in speaking of the 
case in which nuclei of isotopic number 1 are used as targets they state: 
“the kinetic energy of the a-particles, plus the energy of mass of the a- 
particle, and of the nucleus to which it becomes attached, should, in general 
be sufficiently great to incite the emission of a neutron.” 

They predicted, for example, that silicon 29 would react in the following 
way provided an a particle were captured: 


1:50 + 9Het > 4S — 12S + gn! (9) 
since a stable atom of sulphur 32 could be formed by the disintegration of 


the intermediate nucleus sulphur 33. 
They also considered that if the intermediate nucleus formed were phos- 
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phorus 31, which could itself be in a stable state, it would emit a proton 
since it could thus emit its excess energy and return to a stable state. Thus 


Al? + Sa > gP3 — Si + “1H. (10) 


Both of the above have been found to be correct. There is, however, 
another reaction by means of which the excess energy of phosphorus 31 
may be released by giving silicon 30 asa product. This is 


Pa — 13P™ a gn} (11) 
excited unstable 
stateof atom 
stable 
atom 


1-P™ —> 145i + ~ fe (12) 


The use of a star to indicate the energy of activation of the intermediate 
product was introduced by Harkins and Gans® who gave the general reac- 
tion as 

A+B—AB*~C+D (13) 
which may be written 
T+P-I*-C+D (14) 


where a target nucleus (7) struck by a projectile (P) gives an excited inter- 
mediate nucleus J*. In the earlier paper it had already been assumed that 
C or D may be emitted in an excited state, or 


T+P7I*>cC*+D (15) 


which means either the emission of a y-ray or a further disintegration of C*. 

Our use of the star should not be confused with its later use by others, 
who have used it to designate merely the possibility of emission of electrons 
or y-rays, but not the energy of excitation of the intermediate state. 

What has been presented thus far represents the state of the Intermediate 
Nucleus theory at the end of 1935 except that in 1933 the writer determined 
experimentally for reaction (9) that the life of this particular intermediate 
nucleus, with the energies used, is less than 10~’ seconds. This was not 
unexpected as from general considerations it was expected that the life of 
any intermediate nucleus to be less than 10~! seconds, but in this case 
more of the order of 10~'* or 10—"” seconds. 

Nuclear quantum theorists assumed the time of interaction to be of the 
order of R/v, where R is the radius of the nucleus and 1 is the velocity of the 
projectile. However, the theory of the intermediate nucleus assumed that 
the time of disintegration does not depend upon the velocity of the projec- 
tile, except in a very much less direct way. The time of disintegration does, 
however, have a certain indirect dependence upon the energy of the pro- 





r- 
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jectile, as considered later. The theorists also assumed the interaction 
between the projectile and the particles in the nucleus to be much greater 
than that with the whole nucleus, while the intermediate theory assumed 
exactly the opposite. 

Thus the writer considered that both the kinetic energy of the projectile, 
and also the energy of any excess mass which it and the target nucleus 
may possess, are very rapidly distributed in the new nucleus by the motion 
of its component particles, and thus put the nucleus in some excited state, 
the nature of which, together with the possible modes of disintegration, 
decides the life of the particle. 

A general outline of this theory was presented’ in an address given 
January 1, 1936, and the theory was given the support of Bohr® in an ad- 
dress presented January 27, 1936. 

The purpose of this paper is to outline more definitely the postulates of 
the theory. 





PROJECTILE TARGET 


FIGURE 3 
Illustrates the Theory of the Intermediate Nucleus. 


(1) The fundamental theory is that the intermediate nucleus has an 
actual life and a separate existence. 

(2). The life period of any definite intermediate nucleus is independent 
of the components from which it has been formed provided it is in a definite 
excited state, but the state of excitation depends upon both the kinetic 
and mass energies of these components. For example C'**, whether formed 
from (1) B® + H?, (2) B'' + H!' or from (3) Li® + Li®, should have the 
same half-period provided the sum of the mass energy and the kinetic 
energy is the same in all three cases. 

(3) The life period depends upon the particular excited state of J*, and 
on the energy of this state, and upon the nature and the energy states of its 
possible products. 

(4) When any definite intermediate nucleus is formed by the union of 
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two definite nuclei the life period may in general be expected to decrease 
as the energy of the projectile increases. From the quantum viewpoint 
this is due to the general broadening of nuclear levels as the nuclear energy 
increases. On the classical basis it may be explained by the assumption 
that the greater the internal energy of the nucleus, the greater is the proba- 
bility that a single particle [or group of particles such as (mp)2], should in 
any given interval of time, be given enough energy to allow it to escape 
from the nucleus. On the classical basis this would be considered as energy 
of vibration, and the case becomes similar to that of a monomolecular reac- 
tion as considered in chemical kinetics. 

(5) The intermediate nucleus /* may disintegrate into the components 


T and P as well as into other particles. 
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Excitation curve for the action of Helium on Nitrogen. (A) Proton Emis- 
(B) Neutron and Positive Electron Emission (work of Haxel). 
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sion. 

(6) On account of the low decimal mass of helium (0.0034 or 3.2 MEV) 
the emission of helium is favored when the energy of excitation of J* is low, 
the nuclear charge is low, and also a state of even number is involved (Bose 
statistics), since the angular momentum of the a-particle is zero. For pro- 
tons and neutrons the level should be odd (Fermi statistics). 

(7) According to the theory the total number of nuclear disintegrations 
obtained must be equal to the number of syntheses of the excited interme- 
diate nucleus minus the number of these nuclei which return to a stable 
state by y-ray emission alone. Here disintegration into either heavy par- 
ticles, or into negative and positive electrons is considered, but with light 
nuclei it is well known that heavy particles, such as neutrons, protons, a- 
particles, etc., are almost always emitted by the first or primary interme- 


diate nucleus. 
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The only definitely determined case in which an electron is emitted but 
not a heavy particle is the following: 


C2 -+- ~im al ~1N13* (16) 
~ N+ > ~INB* + (17) 
“Ie — IC" 4 “te (18) 


Here, however, the primary excited intermediate product emits a part of 
its energy of excitation and changes into an unstable form. This, which 
has a half-life of 660 seconds, emits a positive electron. 

(8) From (5) the total number of disintegrations is thus in a secondary 
sense dependent upon the nature of the constituents T and P, including 
their energies, but is independent of the nature of the products, C and D, 
formed by the disintegration. 

(9) The types of disintegration which occur and their corresponding 
widths of level depend upon the state of the intermediate nucleus, particu- 
larly upon its energy, and the angular momentum. 

(10) The disintegration of the intermediate nucleus 7* may occur in 
steps, such that a second intermediate nucleus J’* is formed. 

Consider reactions (4), (5) and (5’) for which Harkins and Shadduck 
assumed the existence of an intermediate fluorine nucleus. It is known 
that oxygen 17 is formed by reactions in which (A) hydrogen is formed, 
and (B) a neutron is formed. Now, if the total number of captures of 
helium by nitrogen were known (N,) for any given velocity of the projec- 
tile (helium), this would give the number of disintegrations to give hydro- 
gen (N,,) plus those to give neutrons (N,), 


or N, = N, + N,, at a given energy. (27) 


Let it be assumed that N, is zero for low energies of the projectile, but 
increases rapidly above a range of 3 cm. (4 MEV) for the helium (a). 
It becomes constant above 4 cm. as the range (or energy) increases. Then 
reaction (2) does not begin until the range is 4.5 cm. or more (energy > 
5.5 MEV). The reason for the necessity of this extra energy for (2) is ap- 
parent in the reactions, since in (5) and (5’) a positive electron is formed, 
in addition to the neutron which has a slightly higher mass than the hy- 
drogen of (4). According to (27) as the number of reactions (B) increases, 
the number of (A) should decrease. Now this is just what has been 
found by Haxel (Fig. 4). Thus the study of excitation potentials should 
give evidence for the theory of the intermediate nucleus. Unfortunately 
the total curve for the number of captures is not in general of this simple 
form. 

Let it be assumed that heavy hydrogen is used to synthesize carbon 12, 
by the reaction 








128 CHEMISTRY: W. D. HARKINS Proc. N. A. S. 


5B” + JF C* (19) 
15.0 + (13.7 + 1) = 3.4 + 26.1* MEV + 0.2 MEV. 


and that the hydrogen of mass 2 (deuterium) has a kinetic energy of 1 
million electron volts. The atom itself at rest contains an excess mass 
(mass above a whole number) equal to 0.0147 which corresponds to 13.7 
MEV, while boron 10 has an excess of 15 MEV. The sum of these is 29.7 
MEV. Now the excess energy of a carbon 12 atom in its normal nuclear 
state is 3.4 MEV, so the excited level of C* has an excitation energy of 


O 


PROJECTILE TARGET INTERMEDIATE 
PRODUCT 


/ 


PRODUCTS 
FIGURE 5 
The Disintegration, beginning with the First Intermediate Prod- 
uct may occur in Several Steps or a Disintegration Series. This re- 
quires that the energy of excitation of the intermediate nucleus be 
very high. Such a series is entirely analogous to that exhibited in 
natural radioactivity. 


26.1 million volts, or the star indicates 26.1 MEV, since about 0.2 MEV 
is carried off by the kinetic energy of the carbon 12. 
Now if C'? is formed from light hydrogen 


5B. + ,H!— C* (20) 
11.9 + (7.6 + 10.9) = 3.4 + 26.1* MEV + 1 MEV. 


Here, since about 1 million volts of energy are carried off by the C* as k1- 
netic energy, the hydrogen of mass 1 should be projected with a kinetic 
energy of 10.9 million volts to give a carbon 12 nucleus (C*) of the same 
energy as with a 1 million volt deuteron. 

Provided the same excited state is reached the method of disintegration 
of the excited carbon atom should not be dependent upon its method of 
formation: Thus 

Ci2* — Bu + HH (21) 
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FIGURE 6 

Reactions of Light Atoms according to Theory of the Intermediate Nucleus, ar- 
ranged according to projectiles as named. The target nucleus is represented by a 
black closed circle, the intermediate nucleii by open circles with a symbol inside. m Thus 
+ or — indicates the emission of an electron, p, of a proton, n, of a neutron and a, 
of an a-particle. A black triangle indicates the heavier of the twoend products. The 
nature of the projectile, as for example for an alpha-particle, is given not only by the 
words but also by the fact that the intermediate nucleus lies two squares to the right 
of the target. In the emission of an a-particle the other product lies two squares to 
the left of the intermediate product. If the intermediate product emits a proton the 
arrow shows a unit displacement in the direction of slope —1, that is, the isotopic 
number is increased by 1 but the atomic number changes by —1. Emission of a nega- 
tive electron lowers n by 2 and increases Z by one, for a positive electron the change is 
the opposite. In the emission of a neutron, n decreases by 1. For the emission of 
any particle the change in the plot is just the opposite of that for its addition as a 
projectile. : 











CHEMISTRY: W. D. HARKINS Proc. N. A. S. 


Ce* — Bw 4 FP (22) 
CoP > Be? 4- He* (23) 
C2* > Cu* 4 yl (24) 


are among the reactions which should ensue. In addition 


Cut —> Bu + e+ (25) 
and Be’* — He + He. (26) 


The half-lives of the excited states may be a few seconds, minutes, hours 
or possibly even years if an electron is emitted from an intermediate nu- 
cleus. For the emission of a neutron, proton, helium nucleus or so-called 
heavy particle, the life may be predicted in some cases from the energy of 
excitation of the state. 

Thus Be®* has a narrow energy level from which y-rays alone are emitted, 
and a superimposed broad level from which helium nuclei are emitted. 

From the present value of the energy of this latter level Bethe (private 
communication) estimates the life of this particular intermediate nucleus 
in this state as 10~'* seconds, which is a long life for the emission of a 
“heavy”’ particle from a light nucleus. 

With heavy nuclei and neutrons, energy levels as narrow as 1/50 volt 
have been found in silver. An energy level of this width in an intermediate 
nucleus would correspond also to the general order of 10~'* seconds. With 
a-particles, for example with aluminum, levels of widths of 0.5 MEV are 
not uncommon. 

From this it may be estimated that for the emission of ‘“‘heavy”’ particles 
(particles which are not electrons or -rays) the life periods of intermediate 
nuclei will in general lie between the limits of 10~'* and 10—'® seconds, 
though 10—* and 10~—* seconds are not improbable in certain cases. 

It is of course possible that the life of a nucleus may in some instances 
become so small that the complete independence of the disintegration from 
the method of formation may partly vanish. 

The theory of the intermediate nucleus is illustrated by figure 5, which 
shows the course of most of the known nuclear reactions of light atoms, on 
the basis of this idea. Direct evidence for the theory of the intermediate 
nucleus would be that which indicates that the assumed intermediate, as 
C'**, actually disintegrates in two (or more) steps. In this special case in 
reaction (23) the first helium atom could be emitted with a variety of 
velocities, but in the disintegration of the second intermediate, as Be®*, 
the velocities of the two helium nuclei produced (reaction (26)) must be 
identical when referred to the moving center of gravity of the system of the 
two helium atoms as the origin of the system of coérdinates. Dee and Gil- 
bert® seem to have shown that this is true. 
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